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Abstract 

We study the off-shell mixing and renormalization of flavor-diagonal dimension-5 T- 
and P-odd operators involving quarks, gluons, and photons, including quark electric 
dipole and chromo-electric dipole operators. We present the renormalization matrix to 
one-loop in the MS scheme. We also provide a definition of the quark chromo-electric 
dipole operator in a regularization-independent momentum-subtraction scheme suitable 
for non-perturbative lattice calculations and present the matching coefficients with the 
MS scheme to one-loop in perturbation theory, using both the naive dimensional regu¬ 
larization and ’t Hooft-Veltman prescriptions for 75. 



1 Introduction 


Permanent electric dipole moments (EDMs) of non-degenerate systems violate invariance un¬ 
der parity (P) and time reversal (T), or, equivalently (l], CP, the combination of charge conju¬ 
gation and parity. Given the smallness of Standard Model (SM) CP-violating (CPV) contribu¬ 
tions induced by quark mixing |2| (for a review see [3|), nucleon, nuclear, and atomic/molecular 
EDMs j4j-|s] are very deep probes of the SM 0-term (already constrained at the level of 
9 ~ 1CP 10 ) and of possible new sources of CP violation beyond the Standard Model (BSM). 

In fact, EDMs at the sensitivity level of ongoing and planned experiments probe BSM 
CPV interactions originating at the TeV scale or above (up to hundreds of TeV depending on 
assumptions about the BSM scenario). These new CPV interactions may be a key ingredient 
of relatively low-scale baryogenesis mechanisms such as electroweak baryogenesis (see |9j and 
references therein), making the study of EDMs all the more interesting. EDMs of the nucleon, 
nuclei, and atoms are sensitive to a number of new sources of CP violation, in a complementary 
way [l0|, so that a broad experimental program to search for EDMs in various systems is called 
for (a summary of current status and prospects can be found in Refs. jTTJ[l21). 

Extracting robust information on the new CPV sources from the (non)observation of EDMs 
is a challenging theoretical problem, that involves physics at scales ranging from the TeV (or 
higher) down to the hadronic, nuclear, and atomic scales, depending on the system under 
consideration. The relevant physics at the hadronic and nuclear scale involves strong inter¬ 
actions, and requires the calculation of non-perturbative matrix elements. While interesting 
model-independent statements can be made within a nucleon-level chiral effective theory ap¬ 
proach 13 18 , ultimately the computation of a number of hadronic matrix elements is nec¬ 


essary. Existing calculations of the impact of BSM operators on hadronic EDMs typically 
rely on modeling the strong dynamics in ways consistent with the Quantum Chromodynamics 


(QCD) symmetries, using methods such as QCD sum rules [19-22 and the Dyson-Schwinger 
equations 23,24j (see Refs. |3,[l2|| for reviews). Since models do not rely on systematic ap¬ 


proximations to the strong dynamics of quarks and gluons in the nucleons, current results 
represent in some cases only crude estimates, with different model-calculations differing by 
up to an order of magnitude, depending on the operator under study [l2|. Needless to say, 
this state of affairs greatly dilutes the impact of EDM experimental searches in probing short- 
distance physics. Moreover, the uncertainties affect the robustness of the phenomenological 
studies relating new sources of CP violation to baryogenesis mechanisms (depending on what 
is the dominant mechanism and operator generating the EDM). 

In this context, lattice QCD calculations offer the opportunity to perform systematically 
improvable calculations of the CPV hadronic dynamics. Historically, lattice QCD efforts have 
mostly focused on the determination of the nucleon EDM induced by the SM 9 term |25]-j33|. 
Only recently there has been interest in studying the impact of the leading CP-odd operators 
on the nucleon EDM 


34-36 and the T-odd pion nucleon couplings 37 


This program, however, comes with several challenges, ranging from controlling the signal- 
to-noise ratio on the lattice to studying operator mixing, and matching suitably renormalized 
lattice operators to the minimally-subtracted operators typically used in phenomenological 
applications. In this paper we focus on dehning UV finite CP-odd operators of dimension 
Eve and lower, using a renormalization scheme suitable for implementation on the lattice, and 
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matching this scheme to the perturbative MS scheme to one-loop. 

The paper is organized as follows. In Section [2j we describe the effective theory frame¬ 
work parameterizing BSM effects at low-energy and identify the leading dimension-5 CPV 
operators. In Section [3j we construct the basis of operators needed to study the renormaliza¬ 
tion of the quark chromo-electric dipole moment (CEDM) operator in an off-shell momentum 
subtraction scheme with lion-exceptional momenta. In Section [4| we present the one-loop 
calculations needed to determine the full mixing matrix to 0{a s ) for the operator basis dis¬ 
cussed in Section [3] In Section [oj we give our results for the matrix of renormalization 
constants in the MS scheme, while in Section [6] we specify the renormalization conditions 
that define a regularization-independent (RI) momentum subtraction scheme and provide the 
0(a s ) matching coefficients to the MS scheme. In Section [t| we discuss the consistency of our 
renormalization conditions with the singlet axial Ward identities. In Section [8] we compare our 
results to recent related work |38| that studies the renormalization of the strangeness changing 
cliromo-magnetic quark operator. We end with our conclusions and outlook in Section [9} 

A number of technical issues are discussed in the Appendices. In Appendix [A] we summarize 
our choice of phase convention used to define the CP transformations. The regularization- 
independent calculation is done using off-shell matrix elements with quarks and gluons as 
external states in a fixed gauge. In Appendix [B] we derive the constraints on the mixing 
with gauge dependent and off-shell operators imposed by BRST symmetry. The Peccei-Quinn 
mechanism and its implications for CPV operators are discussed in Appendix[Cj In Appendix [D] 
we discuss the subtleties that arise in the isospin symmetry limit. Finally, in Appendix [Ej we 
summarize the matching coefficients between the MS and the RI scheme. 


2 Framework 


In this section we describe in some detail the hadronic-scale CPV effective Lagrangian induced 
by BSM physics at the high scale. The identification of the CPV combinations of short-distance 
parameters involves several steps. We start our discussion in Section |2.1| by classifying the 


leading BSM-induced operators that can lead to CPV effects at the quark and gluon level. 


We then discuss in Section A2 the relation between CP and chiral symmetry in presence of 
operators that explicitly break chiral symmetry: the CP symmetry that remains unbroken by 
the vacuum takes the standard form given in Appendix [A] only after performing an appropri¬ 
ate chiral rotation of the fields (“vacuum alignment”) that eliminates pion tadpoles |39 m- 


In Section |2.3| we implement the vacuum alignment in presence of higher-dimensional opera¬ 
tors induced by BSM physics and in Section 2.4| we summarize the vacuum-aligned effective 
Lagrangian including operators up to dimension five. 


2.1 CP violation in the Standard Model and beyond 

Assuming the existence of new physics beyond the Standard Model (BSM) at a scale Absm 
v ew , we can parameterize the BSM effects in terms of local operators of dimension five and 
higher, suppressed by powers of the scale A B sm- The new operators are built out of SM 
fields and respect the SU(3)c x SU(2) W x U(l)y gauge symmetries of the SM. The leading 
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CP-violating operators appear only at dimension six 42,43 . Their renormalization group 


evolution from the new physics scale down to the hadronic scale has been studied in several 
papers, most recently in Refs. |44,45], and the resulting effective chiral Lagrangian at the 
hadronic level has been discussed in Refs. 14 161 (for a review see Ref. fl2|). 


In this work, we are primarily interested in the structure of the effective Lagrangian includ¬ 
ing new sources of CP violation below the weak scale. After integrating out the top quark, the 
Higgs boson, and the W ± and Z gauge bosons, the needed operators are invariant under the 
SU(3)c x C/(1)em gauge group. At a scale g < Mw,z, the effective Lagrangian including the 
leading (i.e., originating at dimension six) flavor-conserving CP-violating effects at the quark- 
and gluon-level can be written as follows]]] 


£ cff = £sM| mi=0 - rrii^Li 'ijiRi - ml'ipRi ip Li - ^^flGG 

- 9 1 * * * V \7 e ( rf «- 7) ^LiCTnuF^i)Ri + d G) * ^RiCT^F^ V’n*) 

^ a bsm v ' 

- g(d[ 9) ^LiCr^Ct^ ipRi + d[ 9) * ^Ri<J il vG lLV -fe) 

^ a bsm v ' 

+ + 4-quark operators , (1) 

^BSM 


where e and g are the electric and color charges, v ew is the Higgs VEV (vacuum expectation 
value), the index i runs over the active quark flavors (at /i ~ 1 GeV one has i G {w, d, s}), and 
G^ b = E^ a PG b a p/2. 

The first line in Eq. ([I]) contains the Standard Model dimension-4 operators, including the 
mass matrix put in the standard diagonal form and a common phase, and the QCD 0-term. 
Because of the anomalous Ward identity, a choice of fermion phases can be used to rotate the 
0-term into a CP-odd pseudoscalar quark mass term, instead. 

The second and third lines in Eq. ([Tj) contain the BSM contribution due to the quark 
magnetic (MDM) and electric dipole moment (EDM), and chromo-magnetic (CMDM) and 
chromo-electric dipole moment (CEDM) operators, respectively. Below the weak scale these 
operators are of mass-dimension five: their origin as dimension-6 operators at the high scale is 
hidden in the overall dimensionless factor of n ew /ABSM- It is important to note that the phys¬ 
ical meaning of these operators as CP-violating electric or CP-conserving magnetic moments 
relies on an implicit chiral phase convention. Similar to the quark mass terms, however, these 
operators explicitly break chiral symmetry, thus contributing to vacuum alignment |39|. As 


discussed in Section 2.2, the vacuum alignment, in turn, determines the unbroken CP symme¬ 
try, and a chiral rotation—which mixes the EDM and MDM, as well as CEDM and CMDM 
operators—may be needed to put the symmetry transformation in the standard form. If no 
complex phases appear in the Lagrangian after such a rotation, no physical CP violation can 


1 Without loss of generality, we have performed a SU (uf)l x SU ( hf)r transformation to put the quark mass 

matrix in diagonal form, with complex masses sharing a common phase p, namely rrii = \mi\e lp . Moreover, 

note that the masses to, and 9 include (i) possible threshold corrections, i.e., effects that originate from higher 

dimensional operators, such as H^HGG and iH Hq^q^H , and (ii) corrections induced by mixing with the 

chromo-electric dipole moment at finite quark mass. 
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arise. 


In Section 


2.3 


we, therefore, discuss the combinations of and d^ that are 

independent of such phase choices, and give CP-violating contributions to observables. 

Finally, the fourth line in Eq. ([TJ) contains the CP-odd BSM operators that are genuinely of 
mass-dimension six at low energy, such as the Weinberg three-gluon operator and four-quark 
operators. 

In order to convert experimental results on nucleon and nuclear EDMs into bounds or 
ranges for the short-distance CP-odd couplings, one needs to compute the effect of the CP-odd 
operators in Eq. ([TJ) on hadronic observables, such as the nucleon EDM and the T-odd irNN 
couplings. One essential step in connecting the short-distance physics to hadronic observables 
involves defining UV finite operators in a suitable scheme, whose matrix elements can then 
be computed non-perturbatively using lattice QCD. In this work, we focus on the ultraviolet 
divergences and mixing structure of the leading gauge-invariant CP-odd dimension-5 operators, 
namely the quark CEDM and EDM. These operators are of great phenomenological interest, 
being the leading sources of flavor-diagonal CP violation in several extensions of the SM |3,12|. 
Moreover, since dimension-5 operators can mix only with operators of dimension up to five 
(mixing with lower dimensional operators occurs in mass-dependent renormalization schemes), 
we can consistently ignore operators of dimension six and higher, which we leave for future 
work. 


2.2 CP symmetry and chiral symmetry breaking 

In this subsection we discuss the connection between CP and chiral symmetries. The main 
point is that explicit chiral symmetry breaking selects the vacuum of the theory |39], as well 
as the unbroken CP symmetry. The unbroken CP symmetry takes the standard form given in 
Appendix [A] only after a chiral rotation that eliminates pion tadpoles, he., after implementing 
vacuum alignment |39| discussed in Section 

The CP transformation interchanges left-chiral particles with right-chiral anti-particles. It 
is implemented on chiral fermion fields by 

CV-^lCV = i'y 2 iPl T 

CV-\p R CP = (2) 

where CV is the CP operator (see Appendix [A] for details). The CP operation does not 
commute with chiral rotations, so we can consider its outer automorphisms. In fact, defining 
the chiral rotation operator \ via (i labels quark flavors) 

X X ^L,iX = e _ * Xi/2, 0L j i 

X~ X ^R,iX = e lXi/2 ^R,i , (3) 

one finds 

CV-^L, CV X = ie lXi X2^L T i 

CV- 1 ^ CV X = (4) 


2.3 


4 




where CP X = X~ l CP\. If chiral symmetry is a good symmetry of the Lagrangian £ 0 , then 
each of these is an equivalent CP symmetry. 

Because of the spontaneous breaking of chiral symmetry, almost all the CV X are sponta¬ 
neously broken by the vacuum of the theory. In this case, it is convenient to make a chiral 
phase choice such that the vacuum has a zero expectation value for all the flavor bilinears of 
the form (^ 75 ^). In fact, it is only with this phase choice that the pions, the Goldstone 
modes of the broken chiral symmetry, correspond to the operator ^ 75 ^. With this choice of 
phases, in the “reference vacuum”, the CP symmetry CP 0 stays unbroken by the vacuum; we 
implicitly make this choice throughout this paper. 

We next consider the effect of explicit chiral symmetry breaking. For a small explicit 
breaking of chiral symmetry, encoded in a new term 5C in the Lagrangian £ = £ 0 + £ 5C 
(with e <C 1 ), chiral perturbation theory is expected to be a good guide to understanding 
the structure of the theory. But, because of the explicit breaking of the chiral symmetry, the 
vacuum is no longer degenerate: the explicit breaking chooses a direction in chiral space with 
which the vacuum aligns 139]. If this does not match the “reference vacuum”, large corrections 
appear due to degenerate perturbation theory. 

To avoid this problem, it is convenient to perform a chiral transformation y so that the 
explicit chiral symmetry breaking 5C selects the reference vacuum, in which the unbroken CP 
symmetry takes the standard form, namely CP 0 . The way to do this is to impose the condition 
that the vacuum state does not mix with the Goldstone state 


39-41 , i.e., 


<7r|<$£|fi> = 0, 


(5) 


where <5£ are the chiral breaking terms after such a rotation and Q) and | 7 r) are the reference 
vacuum and Goldstone pion states respectively. If the only chiral breaking comes from the 
mass terms, this can be accomplished by rotating away the flavor non-singlet CP-violating 
mass terms in Eq. (|T|) by the appropriate chiral transformation y. 


2.3 Vacuum alignment in presence of higher-dimensional operators 

We now discuss the vacuum alignment in presence of higher-dimensional operators induced 
by BSM physics. After a general discussion of the chiral transformation needed to enforce 
Eq. (J5|, we specialize to the case in which the dominant source of chiral symmetry breaking is 
provided by the quark masses, and the dominant BSM operators are the quark (C)EDM and 
(C)MDM. In this case we present the vacuum-aligned effective Lagrangian in both scenarios 
with and without the Peccei-Quinn (PQ) mechanism |46|. 

Except for GG, all terms of dimension five and lower in the Lagrangian defined in Eq. (Jl| 
that violate CP are fermion bilinears that also violate chiral symmetry. Each of these terms 
mixes with a CP-conserving one under chiral rotation, and it is conventional to treat the two 
as real and imaginary parts of a single operator. Generalizing Eq. (|TJ) let us write the chiral 
and CP-violating part of the Lagrangian involving quark bilinears as 

5£ = - [df Of + h.c.] = - [ Re df Re Of + Im df Im Of] ( 6 ) 

i,a i,a 
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where 


Of = 


Re Of — Of + O 


af 


Im Of = i 


Of - O 


af 


( 7 ) 


i is a flavor index and a parameterizes the different operators, characterized by the structure 
r°. The first few operators are the mass term (a = 0), the quark CEDM (a = 1) and the 
quark EDM (a = 2) 


+ h.c. = ilJi[(Remi) + 

7l ^ • ^ ~ - Vew g ^[(Red ( ^ 

l BSM 


d]0] + h.c. = g A (R .ed\ 9> )a flu G flu + i(Imd ( : 9) )a lxu G lil/ 'y 5 


d^Oi + h.c. = 


v. 


2A 2 

zii BSM 


e A 


(Red^Wi^ + z(Imd™)^^75 


(7) > 


A 

A 


In this notation, under a chiral rotation x (parameterized by Xi 


( 8 ) 

(9) 

( 10 ) 


df -2 dfe iXi 

@ @ + X l + ' ' ' + Xn F , ( 11 ) 

and we seek a chiral rotation such that Eq. ([5]) holds and at the same time 0 —> 0. 

To implement Eq. ([5]), we need to introduce the non-perturbative matrix elements 


Afj = Aj | Im Of\Q) ( 12 ) 

where the state 717 ) is interpolated by the field Then the mixing of the vacuum with 

the neutral Goldstone modes (17Tj) — \nk))/V^ is proportional to Y2i a Im df (A— A f k ). The 
condition in Eq. (|5]) for each neutral Goldstone mode (17r ^ ) — |7Tfc))/ a/2 becomes 

E Im «e ix -) [A“ - A“J = 0 , k = l,j = 2,...,n F , (13) 


Since the unperturbed Lagrangian Cq is SU(uf)v symmetric, the matrix elements can be 
written in terms of two constants, the diagonal Af and the off-diagonal A^, defined by A“- = 

?(i 


A fjSij + A 


Sij). Eq. (13) implies, for each flavor i = 1 


Im (dfe 1Xi ) 


r ^ = k 


(14) 


where r^ a> = (Ag — Ay)/(A(j — Ay) (we divided out the matrix elements of the dimension-3 
operator Im 0 { - = mid k is a flavor-independent constant. Defining 


di = \di\e** = Y2 d > (a \ 

a 

the chiral rotation we want needs to satisfy, for each i, 


(15) 


|dj|sm(Xi + ^i) = k ■ 


(16) 
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Moreover, to implement 9 —y 0, one needs 0 + Y2iXi — 0? or, equivalently, the constant k needs 
to satisfy 


0 - V’ + E. sin x ) — 0. 


(17) 


i 


i 


Eqs. (16) and ( [l7| ) provide a system of equations for Xi and which does not have a closed 
form solution for n E > 2. On making the chiral transformation dictated by Eqs. (16) and © 
we find that CP violation is proportional to 


SC 


cpv — 


E 

i,a 

E 


k Re + yJ\di\ 2 -K 2 Im % 


d. 


di 


Im Of 


-d{6 - 0tot) Re C -j- + \di \Im 

di d{ 


Im Of, 


(18) 


where d 1 = |dj| 1 and 0 t ot = 'Yhi ( t>ii an d the second line is obtained by solving Eq. © 
for small n/\di\, which is appropriate when 9 is small and the dominant chiral violation comes 
from a real mass term (the latter condition implies 0* 1). 

Notice that if there is a single operator Of that is the only source of CP violation, then 
di oc df, and the second term is zero. This is because in this case this term is also the only term 
that explicitly breaks the chiral symmetry and the vacuum aligns itself with this direction. 
As a result, performing a chiral rotation to make the vacuum have the conventional chiral 
phase removes any imaginary part from the operator, and CP violation can only come from 
the anomalous chiral rotations. In this case, however, the CP violation is proportional to the 
harmonic sum of the chiral violations from each flavor, and therefore vanishes if any flavor 
remains chirally symmetric. 

In what follows, we will instead consider the situation where the dominant chiral breaking 
is always due to the a = 0 mass term, i.e., di oc approximately, and consider the case where 
all flavors are massive. Only in this case, the dominant source of CP violation is proportional 
to Im df. Consistent with this assumption, when studying mixing and renormalization we will 
keep in the operator basis terms proportional to the quark mass matrix. 

With these assumptions and after vacuum alignment, the explicit form of Eq. (18), spe¬ 
cialized to the case of a Lagrangian containing a mass term, quark EDM, and quark CEDM, 
is 


5C C pv = {o - ^Tr [M 1 ([dcs] - mJM 1 [^cm])]) 

+ 2^75 ([ d C E] - {d CM ]) 

- y >75 ([ d CE \ - mjM~ l [ d C M ]) 

qp _ _ 

- j F{{d E } - m+QM.- 1 [d M ]) ^ , (19) 

where we defined r = -A 1 ) = (A<j — Ay)/(A| — Ay) and neglected A 2 ) = 0 (ckem r^ 1 )). We 
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further defined 



lu\ 


^ m u 0 0 ^ 

i) = 

d 

M = 

0 m d 0 


W 


V 0 0 m s j 


and the matrix-valued CEDM and CMDM couplings as 


(20) 


[dcE] 


^ew 

A2 

iv BSM 


(Im du ^ 0 

0 Im 

\ 0 0 


0 > 

0 

Im di 9) ) 


[' dcM\ 


hw 

T 2 

iv BSM 


/ Re d i9) 0 0 \ 

0 Re d!f 0 
^ 0 0 Re di 9 '*) 


1 


( 21 ) 

with analogous dehnitions for the electric [d#] and magnetic [d^] couplings. Finally, 6 = 
6 — ufP with ufP the phase of the determinant of the mass matrix before the anomalous 
chiral rotation renders it real, and m* is the reduced quark mass 


m* = 


m s m d m u 


m s (m u + m d ) + m u m d 


( 22 ) 


The first term in Eq. (19) is the familiar 6 term, shifted by a correction proportional to the 


quark CEDM and a second correction, proportional to the coefficients of the CMDM multiplied 


by 9. The third and fourth lines of Eq. (19) contains the quark (C)EDM operators, which after 


vacuum alignment receive a correction proportional to the (C)MDM coefficient multiplied by 
6. Moreover, vacuum alignment causes the appearance of a complex mass term, proportional 


to the same combination of the CEDM and CMDM coefficients (second line of Eq. (19)). 

The above discussion is valid in absence of PQ mechanism |46|. As we review in Ap¬ 
pendix [Cj if CP violation arises only from the mass term, the PQ mechanism dynamically 
relaxes 6 to zero. In the presence of other CP-violating sources, like the quark CEDM, the 
Peccei-Quinn (PQ) mechanism causes 9 to relax to a non-zero value d in d, proportional to the 
new source of CP violation. In particular, as we discuss in further detail in Appendix [CJ in 
the presence of the quark CEDM 


«i»d = j'Tr [AT 1 [d CE }} , 


(23) 


thus enforcing a cancellation between the first two terms in Eq. (19). Since #i n d is suppressed 


by two powers of A B sm, terms proportional to 9[dcM\ hi Eq. (19) become effectively dimension 
eight, and can be neglected. Thus, if the PQ mechanism is at work, the first line of Eq. (19) 
vanishes and the terms proportional to 6 in the second and third line of Eq. (19) can be 
neglected, leading to 


5Ccp V = ^ ^75 [d C E\ ^ ~ y ^^>75 G^ [< d C E\ if* - y 75 F 9U [dp] if* , (24) 

with both CEDM and pseudoscalar quark density with flavor structure dictated by \dcE\- 


Eq. (19) and Eq. (24) provide the vacuum-aligned low-energy Lagrangians, in presence of 


BSM sources of CP and chiral symmetry violation. They are particularly useful within the 
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chiral perturbation theory framework, as they guarantee the cancellation of tadpole diagrams 
in which Goldstone modes are absorbed by the vacuum. This form of the CP-violating per¬ 
turbation allows one to identify what non-perturbative matrix elements are needed in order 
to address the impact of a BSM-induced CEDM operator on the nucleon EDM, i.e., the de¬ 
pendence of d n on \dcE\- Both with and without PQ mechanism the effective Lagrangian 
involves the CEDM operator as well as flavor singlet and non-singlet pseudoscalar quark op¬ 
erators. Moreover, at the lowest order, the effect of flavor non-singlet operators 

is proportional to insertions of the flavor-singlet density 'ipijs'ijj. This is very simple to see 
within the functional integral approach, in which can be eliminated through a non- 

anomalous axial rotation. The same result can be obtained within an operator approach. In 
this framework, using soft-pion techniques, one can show that a cancellation occurs between 
non-tadpole and tadpole diagrams with insertion of leaving a term proportional to 

the insertion of In absence of PQ mechanism, the resulting flavor-singlet pseudoscalar 


insertion proportional to [cIce] cancels exactly the existing singlet term in Eq. (19). If the PQ 
mechanism is operative, the resulting flavor-singlet pseudoscalar insertion is proportional to 
mj ind . The net effect is equivalent to replacing Eqs. (|I9| and (24) with 


— — XQ — — 

5C C pv = m* 6 - -j- {[dcs] ~ mJM -1 [ d C M ]) V’ 


le 


SjC pq 

0L 'CPV 


= m* Q- md i 75 ^ 


^cr^ 5 F^ iU {{d E } -m*0M 1 [d M ]) (25a) 

^ V’CV 75 G^ v [< d C E} ~ ^ $cr ^75 F^ [d E \ ^ ■ (25b) 


These can be regarded as partially aligned effective Lagrangians, in which only the domi¬ 
nant mass term has been aligned to eliminate pion tadpoles, while the BSM perturbation is 
not aligned. While the physics cannot depend on the choice of equivalent parameterization 


Eq. (19), Eq. (24) and Eq. (25), use of different effective Lagrangians is a matter of con¬ 


venience, depending on the non-perturbative approach employed to study hadronic physics. 


Starting from the Lagrangian in Eq. (25), in the chiral effective theory approach tadpole di 
agrams arise, that can be dealt within perturbation theory fl4 


On the other hand, in a 
non-perturbative approach based on the functional integral, such as lattice QCD, the partially 
aligned Lagrangian can be more convenient: it shows that the only needed non-perturbative 
matrix elements involve the (C)EDM operator and the singlet pseudoscalar density (or equiv¬ 
alently GG). 


2.4 CP-violating effective Lagrangian at the hadronic scale 

To summarize the above discussion, at the hadronic scale (// ~ 1 GeV) the vacuum-aligned 
flavor-conserving effective Lagrangian including the leading BSM sources of CP violation (up 
to dimension five) can be written as follows, 

£ = £qcd+qed - - $ [8M\ 

- ^ i’cr^F^ [De] Q Il> ~ Y [ Dce\ ^ , (26) 
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where 


Q = 


(q u o o ^ 
0 q d 0 
v 0 0 q s J 


(27) 


Here we are neglecting operators that are total derivatives and/or vanish by using the equations 
of motion (EOM), needed later on when we impose off-shell renormalization conditions at finite 
momentum insertion. The matrix-valued CP-violating couplings [hA4], [Dce\, [De\ are related 
to the short-distance couplings of Eq. ([Tj) via Eq. (21) and Eq. (19) or Eq. (24), depending on 
whether or not the PQ mechanism is assumed. The pseudoscalar mass term [AM] in general 
has a non-singlet structure in flavor space, though at leading order, its physical effects can be 


related to a flavor-singlet mass term as discussed in Section 2.3 (see Eqs. (25)). 


3 CP-odd operators of dimension < 5 

The only T-odd and P-odd operators of dimension Eve appearing in the low-energy effective 


Lagrangian Eq. (26) are the quark EDM and CEDM, whose mixing and renormalization we 
wish to discuss. 

The analysis of the quark EDM is relatively simple: this operator is a quark bilinear 
from the point of view of strong interactions, and it is simply related to the tensor density. 
Knowledge of the nucleon tensor charges immediately allows one to extract the contribution 
of the quark EDM to the nucleon EDM |47, 48 . To lowest (zeroth) order in electroweak 


interactions, this operator renormalizes diagonally, precisely as the tensor density. Since we 
are not interested in the hadronic matrix elements to a precision of order olem/w < 1%, we 
neglect the quark EDM mixing with any other operator. 

On the other hand, the quark CEDM operator does not renormalize diagonally: it mixes 
with the quark EDM and other operators of dimension five or lower. The mixing structure is 
particularly rich if one considers renormalization within a so-called regularization-independent 
(RI), momentum subtraction (MOM) scheme, amenable to non-perturbative calculations in 
lattice QCD |49|. In this family of schemes the renormalization conditions are imposed on 
off-shell quark matrix elements in a fixed gauge, thus requiring the inclusion of operators that 
do not contribute to physical matrix elements, such as total derivatives and operators that 
vanish on-shell by using the equations of motion (EOM). We next discuss the relevant operator 
basis, the mixing structure, and the strategy to determine the renormalization matrix. 


3.1 Operator basis 

The implementation of RI momentum subtraction schemes requires working in a fixed gauge. 
With gauge fixing, full gauge invariance is lost and the action is only invariant under BRST 
transformations 50,5ll. A given gauge invariant operator O (we have in mind the quark 


CEDM) mixes under renormalization with two classes of operators of the same (or lower) 
dimension 1 52,531 : (i) gauge-invariant and ghost-free operators with the same symmetry prop¬ 
erties as O (Lorentz, CP, P) that do not vanish by the EOM; (ii) “nuisance” operators allowed 
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by the solution to the Ward Identities associated with the BRST symmetry: these vanish by 
the EOM and need not be gauge invariant. The “nuisance” operators can be constructed as 
off-shell BRST variation of operators that have ghost number —1, but otherwise with same 


symmetry properties as O, as discussed in Ref. 52 and detailed in Appendix E J 


Following the above general prescription, we have constructed the basis of CP-odd (T- and 
P-odd) operators that mix with the quark CEDM operator (the CP transformation properties 
of fields are reviewed in Appendix 0- We present our results for rip = 3. To restrict the 
possible structures in flavor space we use the spurion method. While the effective Lagrangian 
in Eq. (26) is not invariant under chiral transformations on the quark fields i^l,r —> Ul,r^l,r 
with Ul,r G SU(3)l,r, one can formally recover chiral invariance by assigning spurion trans¬ 
formation properties to the CEDM coupling matrix {[Dce] —^ Ul[Dce]U ^), the mass matrix 
(Ad — >■ UlMUr), and the charge matrix (Q — * Uj J jiQU\ j R ). One then includes in the basis 
operators that are chirally invariant in the spurion sense, and are linear in the CEDM spurion 
[D C e]- Eventually, we set [Dce] —> t a {a = 0,3,8), where t° = l/\/6I 3x3 is proportional to 
the identity matrix in flavor space, while for a = 3, 8 , t a = A a / 2 , with A a the SU(3) Gell-Mann 
matrices (normalizations are such that Ti 'p(t a t a ) = 1/2 for a = 0,3,8). 

In our basis we include operators proportional to the quark mass matrix for two reasons: 
(i) the identification of the CPV terms in Eq. (26) assumed the quark mass to be the dominant 
source of explicit chiral symmetry breaking; (ii) we wish to include the effect of the strange 
quark, for which m s jA qcd is not a big suppression parameter. 

Finally, in order to present the operators that vanish by the EOM in a compact form, we 
introduce the combinations: 


fpE = (iD» 7 m - M)i/j , 

D„ = d„~ igA‘X - ieQA W 

(28) 

4>e = + M) , 

= K + vA a „ T ‘ + “CA’ • 

(29) 


Note that i^e transforms under CP in the same way as %jj (see Appendix |A| . 

Next, we enumerate the operators of dimension Eve and lower that can mix with the quark 
CEDM: 

C = ig , (30) 

labeled by the flavor-diagonal structure t° (a = 0,3,8). We will use the notation O,- to 
indicate the i th operator of dimension d. If the regularization breaks chiral symmetry, i.e. an 
additional left-right spurion (proportional to the identity in the case of Wilson fermions) is 
present in the effective Lagrangian, the CEDM operator can mix with additional operators. 
While we will restrict our analysis to the case of good chiral symmetry (which can be attained 
on the lattice by using domain-wall (55] or overlap (56 fermions), we will nonetheless identify 
the additional operators appearing at a given dimensionality. Finally, note that there are no 
CP-odd operators containing ghost-antighost fields up to and including dimension five. 

2 See Ref. j54) for an application of this formalism to the CP-even sector of QCD. There is a one-to-one 
correspondence between our operator basis and the one of Ref. [54], provided we drop the total-derivative 
operators from our basis and set m = 0, as done in Ref. (541. 
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3.1.1 Dimension 3 

At dimension three there is only one operator allowed by the symmetries: 

0 (3) = P = 0i 75 t a 0 . (31) 

This operator mixes with the quark CEDM even in the absence of other sources of chiral 
symmetry breaking, such as mass terms or regularization artifacts. Therefore, the lattice 
operator Cl requires subtraction of power divergences due to mixing with the lower dimen¬ 
sional operator Defining the subtracted operator C = Cl — ZPl, one can determine Z 
by requiring that the quark two-point function vanishes at a given symmetric kinematic point 

p 2 _ p /2 _ ^2 _ _^2 £ or rn ^ —> 0, namely Tr (vf^t 0, ) = 0. 

V / a 0 

3.1.2 Dimension 4 

Assuming good chiral symmetry, there are no dimension-4 operators that mix with the quark 
CEDM operator. If the regularization breaks chiral symmetry in a flavor blind fashion, the 
CEDM can mix with the following operators: 

GG , (-07^75^-0), 0^75 {M,t a } 0, Tr [A4C] 0^750, Tr [A4] 0i7 5 t a 0. (32) 

3.1.3 Dimension 5 

At dimension Eve, fourteen Hermitean operators are present. The first ten operators are 
gauge-invariant and do not vanish by the EOM. The latter four are “nuisance” operators. To 


all operators we assign a number and also a more suggestive name: 

of = C = ig 0 ^G 0 X 0 ^ = \ {o'* 75 + 75^) (33) 

Of = d 2 P = d 2 (0i 75 t“0) (34) 

Of = E = ^ t“}0 (35) 

Of = ( mFF) = Tr [MQ 2 t a ] ^ a0 F lw F a g (36) 

Of = (m GG) = Tr [MP] \e^GfGf (37) 

Of = (md-A )i = Tr [Mt a ] d 0 (07 M 7 5 0) (38) 

Of = (m d ■ A) 2 = ^ c0 (07^75 {M, P} 0) - ^Tr [Mt a ] c0 (07 M 7 5 0) (39) 


3 Since dimensionally regularized operators do not mix with lower dimensional operators at any finite order 
in perturbation theory, we will, when necessary, use a subscript L for operators regularized in a scheme, like 
the lattice, that includes a hard cutoff. 
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f = {m 2 P) 1 = ^ ^*75 {M 2 , t a ] -0 

(40) 

'g 5) = ( m 2 P) 2 = Tr [M 2 ] 

(41) 

l !o' ) = (m 2 P) 3 = Tr [Mt a \ 

(42) 

11^ = Pee = ^Elbt a '^E 

(43) 

'12 = d ■ A e = 9 /i [-0 i? 7 M 75t a '0 + 5t a ^E] 

(44) 

'13 = A d = ^75 -$E$ 75 t a ^ 

(45) 

'u = = Y (ip{Q, P}/ (7) 75-0e - ^e{Q,C}4 (7) 75^) ■ 

(46) 


With a flavor blind breaking of chiral symmetry, the CEDM can mix with additional 
dimension-5 operators, namely: 


Tr [M\ <9 /i ('07 M 75fh/’), Tr [M] , (TibW) 2 'ipi'y^ip, 

Tr [M. 2 t a ] 'ftijs'i/j, Tr [M\ Tr [Mt a ] . 


(47) 


In the perturbative analysis presented below, we will use dimensional regularization. For 75 
we will present results for both the naive anti-commuting scheme known as naive dimensional 
regularization (NDR) and the consistent’t Hooft-Veltman (HV) scheme (see [53 and references 
therein). It is important that the regulator does not break the hermiticity of the operator basis: 
when considering operator insertions in the dimensionally regulated theory, care must be taken 
to ensure that the operators remain Hermitean for arbitrary space-time dimension d. This is 
essential in order to obtain correct results for the finite parts of the diagrams. In what follows 
we will need to insert in loop diagrams, so we provide in Eq. (33) the explicit Hermitean 
form of ()\ , valid both in HV and NDR schemes. 


3.2 Mixing structure and Renormalization Conventions 

The relation between renormalized operators (0*) in any given scheme and bare operators 
(Oj- 0 '*) (expressed in terms of the bare fields) can be written as: 

of 1 = Z« O, . (48) 

The renormalization mixing matrix Z i3 is scheme-dependent and has the general structure 
given in Table [l][^J This structure is dictated by several considerations, including (i) power¬ 
counting (some operators are effectively of dimension three and four with either factors of 
masses or external derivatives and cannot mix with genuinely dimension-5 operators); (ii) 
BRST invariance [521; (iii) vanishing by EOM or at zero four-momentum injection. Indicating 

4 Working to first order in insertions of the new physics operator, each sector labeled by the diagonal flavor 
structure t a (a = 0,3,8) renormalizes independently, so that the renormalization matrix has a block-diagonal 
form in flavor space. 


13 




the gauge-invariant operators that do not vanish on using the EOM (0 ? - 5 ' ) for i — 1,10) by 
O a and the “nuisance” operators (Op' ) for i = 11,..., 14) by N a , the renormalization matrix 
has the block-structure 


o(°) 

AT(°) 


0 




(49) 


The divergent part of (proportional to l/(d — 4) in dimensional regularization or log A 2 in 
a cutoff theory), controlling the physical anomalous dimension, is independent of the gauge¬ 
fixing choice [521. 

In the following we will provide Zo in the MS scheme and in a momentum subtraction 
scheme to one loop order. We will perform the calculations in dimensional regularization 
(d = 4 — 2e) and will present results in both the HV and NDR schemes |53| for 75 and 
the 7-matrix algebra (we use the definition of 75 and e^ ua P given in Ref. m- To extract 
the operator renormalization matrix we define the held, coupling, and mass renormalization 
constants as 


V> (0) = \/z^^ 


(50a) 

Af=s/ziA r 


(50b) 

g {0) = Z a 9 ^us 

g7fi/ 2 


VMS ~ V ( 47r ) 1/2 

(50c) 

?7r 0) = Z m m . 


(50d) 


Here, as usual, /j denotes an arbitrary parameter with dimensions of mass, introduced to 
keep the renormalized coupling g dimensionless ([$] = 0), while [m] = 1, [ip\ = 3/2 — e, and 
[Ay\ = 1 — 6. Note that g and ot s = g 2 /(47r) depend on both /i and e, so that da s /d( log/r) = 
-2ea s + 0(ct 2 ). 

Finally, let us discuss different conventions for the renormalization factors for fields, cou¬ 
plings, masses, and operators, generically denoted by Z. Our definitions in Eqs. (48) and ( |50| ) 
follow the notation typically used in the perturbative QCD literature (see for example |58|). 
However, we warn the reader that the lattice community typically uses a different convention 
(fleshed out explicitly in Ref. [59]), which is related to the one followed here by replacing 
everywhere Z —> Z~ x . 


4 Green’s function calculations 

In order to determine Z VJ and the relation between MS and the RI-SMOM scheme to be defined 
in Section [6] below, we will study amputated two- and three-point function^] with operator 
insertion. These are shown in Fig. [l]and defined as follows: 

J d A x e~ iq ' x (g(p',e*')\0(x) \g(p,e)) = (2vr) 4 5 {A \q + p-p) e*V) T%(p,p') e u {p) (51) 

5 Since the terminology of lattice simulations also counts the points at which the operator is inserted, these 
correspond to three- and four- point functions in that terminology. 
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Table 1: Mixing structure of the dimension-5 operators, with “x” representing non-zero entries. 
Throughout, we neglect effects proportional to the electroweak coupling chew- 


d 4 x e iqx (q(p')\0(x) \q(p)) = (27t) 4 8 ( ' 4 \q + p — p l ) u(p , )Yq\p,p')u(p) 


(52) 


J d 4 xe iq ' x {q(p'), g(k, e*)| 0(x) \q(p)) = (27 k) 4 8^ 4 \q + p — p' — k) w(p')T® (p,p', fc)'ix(p). (53) 

To minimize notational clutter in the above equations and throughout the paper we will 
suppress the color indices, which can be restored as follows. The gluon two-point function 
Y^q carries the color structure 8 CC , where c, d are the octet color indices labeling the two 
amputated gluon external legs. The quark two-point function T)y carries the color structure 
Sij, where i,j are the color indices labeling the two amputated quark external legs. The 
quark-quark-gluon three-point function T^"* carries the color structure T*-, where b is the 
octet color index labeling the amputated gluon external leg and i , j are the color indices 
labeling the amputated quark external legs. Moreover, in our notation Tq is linear in the 
gluon polarization vector, i.e. Y ( q = e b *(k) T^A Analogous definitions exist for the photon 
two-point function and the quark-quark-photon three-point function, which we will denote 


by Y t ^ ( ' 1 \p ) p > ) and T [ Q n \p,p',k) (the latter carries color structure 5^). The assignment of 
momentum flow in these two- and three-point functions with operator insertion is shown in 

Fig. [lj 


-1(3,7) 1 


In any scheme, the renormalization factors Z VJ introduced in Eq. (48) can then be deter- 
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Figure 1: Momentum flow of generic diagrams contributing to the quark-quark, gluon-gluon, 
and quark-quark-gluon Green’s functions with operator insertion. The shaded blob represents 
the operator insertion with incoming 4-momentum q and higher order corrections. In the 
four-point function, the gluon (photon) momentum is labeled by k. 


mined by imposing conditions on the two- and three-point functions defined above. Working 
to first order in cn s , the needed Green’s functions with insertion of 0\ = C read: 


p( 2 ) _ ~p( 2 ) 

1 C - 1 c 


1—loop 




( 2 ) 

(5) 


j T^ 1 


tree 


(54) 


~p(3) _ p(3) 

1 C ~ 1 C 


1—loop 


+ (z^Zg^iz- 1 ) n -i)r 


43) 

c 


tree 


E 




tree 


(55) 


p/zz/ _ p/zi/ 

C C 


1-loop “ u : , 


tree 


(56) 


The simplest perturbative scheme is MS, in which one determines the (Z 1 )i J - by requiring 
cancellation of the poles in e = (4 — d)/2. Similar relations involving insertions of allow 
one to determine the remaining entries Z^ of the renormalization matrix. This program 
requires computing the tree-level and one-loop results for the two- and three-point functions, 
to which we turn next. 


4.1 Tree level matrix elements 

In this section we give tree-level results for the gluon two-point functions Y 1 q (p,p'), the quark 
two-point functions Tq\p,p'), the gluon-quark-quark three-point functions VQ\k,p, p'), and 
the photon-quark-quark three-point functions TQ ,7 \k,p, p'), for all the relevant operators 0\ d \ 
The only operator with non-zero two-gluon matrix element at tree level is 0\' = mGG: 

(p, p') = Tr \Mt‘) x 4 r'p.pj . (57) 

U 5 

An analogous result holds for the photon two-point function Y lll/ ( ^\p, p'). 

o 4 

In Tables [2j |3j and |1| we give the tree-level 1-particle irreducible (1PI) matrix elements 
T o\P’P')> r o\k,p,p'), and Y^Q l \k^p,p t ) for each operator. Throughout, we use the notation: 

a(a,b) = a^a^b u 
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eja,b,c) = . 


(58) 


Finally, for a given operator 0, non-lPI tree level contributions to the three-point functions 
(see Fig. [2J can be expressed in terms of quark and gluon two-point functions as follows 


F o\PiP'i k ) = ~9<f ^ + ^ + 2 m r o } (p, k + p) ~9^o\p- k iP') r T '1 9 "V 


+ m 


s — 


u — rrr 


- f r o(p-v\k)fl , 


(59) 


where s = (p 7 + A;) 2 , u = (p — k ) 2 , t — (p' — p) 2 . 


0 

~p( 2 ) 

1 o 

0 ( 3 ) = P 

i~/ 5 t a 

of 5) = d 2 p 

Og^ = (md ■ A) i 
Of 5) = (md ■ A) 2 
0 ^ 5) = (m 2 P) i 
Oi 5) = ( m 2 P ) 2 
Of? = (m 2 P) 3 

Of? = Pee 
Of? = 9 • 2b 

of? = 

-iq 2 ^t a 

Tr [Mt a ] z^7 5 

(| {M,t a } - |Tr [. Mt a ]) ii 75 

| {At 2 ,r}*7 5 

Tr [At 2 f a ] *75 

Tr [Atf a ] iMjs 

-i p ■ p't a - ^{M 2 ,t a } + icr(p,p')t a + |{At,t a }^ 75 

i q 2 t a — (At, t a } <fr — 2 ia(p,p') t a 75 

—i (p 2 + p ,2 )f a - \ {M,t a } ^ 75 


Table 2: Non-vanishing tree-level 2-point functions with operator insertion. For notational 
conventions and momentum flow, see discussion below Eq. (53). 


4.2 One-loop Green’s functions with CEDM insertion 

At one-loop level, we regulate the diagrams with dimensional regularization, following the 
notation introduced in Section 3.2 Working in general covariant gauge (with gauge fixing 
parameter £) [ffj we have computed both the divergent and finite parts of the Green’s functions 

6 Feynman gauge corresponds to £ = 1, while Landau gauge corresponds to £ = 0. 
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0 

rg 3) (ipi) 

of = c 

Of — Pee 
of = d-A E 

of = A d 

2ga(e*,k)^ 5 t a 

—ig e* ■ (p + p') — ia(e*,p — p') 7 5 t a 

2 g cr(e*,q)^ 5 t a 

—ig e* ■ {p + p') + ia(e* ,p — p' — 2k) 7 5 t a 


Table 3: Non-vanishing tree-level 1PI quark-quark-gluon 3-point functions. For notational 
conventions and momentum flow, see discussion below Eq. (53). 


0 

rg ,7) (ipi) 

Of = E 



e{Q,t a } cr(e*,k )75 



Of = Pee 


e* ■ (p + p') — ia(e*,p — p') 

75 


of = d-A E 



e{Q,t a } v(e*,q) 7 s 



Of = A d 

-f {QX} 

e* 

• {p + p') + ia(e* ,p — p' — 2k) 

75 

P)( 5 ) — A 

C^14 — r) 

-f {QX} 

e* ■ (p + p') — ia(e*,p — p') 

75 









Table 4: Non-vanishing tree-level 1PI quark-quark-photon 3-point functions. For notational 
conventions and momentum flow, see discussion below Eq. (53). 



Figure 2: Non-lPI diagrams contributing to the quark three-point function. The shaded blob 
represents the 1PI contribution to the relevant two-point function with operator insertion. 



Figure 3: Diagrams contributing to the quark two-point function. The dot denotes the inser¬ 
tion of the CEDM operator. 
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at generic kinematic points, before specializing to non-exceptional momentum configurations 
needed to define the operators in the RI-SMOM scheme (see Section [6]). Specifically, for the 
two-point functions (with p+q = p') we work at the symmetric point p 2 = p' 2 = q 2 = —A 2 . For 
the three-point functions (p + q = p' + k) we work at the non-symmetric point S characterized 
by p 2 = p' 2 = k 2 = q 2 = s = u = i/2 = —A 2 . We will provide the motivation behind this 
choice in Section [6] 

Throughout this work we will denote the SU(Nq ) color factors as follows: 

TV 2 — 1 1 

° F = ~2nT' Ca = Nc ’ Tf= 2’ (60) 


4.2.1 Quark two-point function 

At one loop, T c\p,p') receives contributions from the diagrams in Fig. [i]and reads: 


r c\p,p) = ^{ (p 2 + P 2 ) l5t a 


xa 


An r 

+ 75 

+ 75 


3 Cf ( - + log ) + /o 


ZC F (1 

— ( J + ‘° g A5 +/‘+°( a? 


P 


-sc F [- e +io g ^\+ f 2 + °r^ 


m: 


(61) 


where 


cHV 

Jo 

22 

fNDR ^ 

;° - 3 

x 3 C F 

(62) 

f HV 

J 1 

= —3 C F , 

fN D R = _2 

J1 3 

3 C F 

(63) 

rHV 

J 2 

= -y x 3 C F 

rNDR 

? / 2 — 

— 3 x 3 C F . 

(64) 


4.2.2 Gluon two-point function 


As illustrated in Fig. [4| at one loop level three diagrams contribute to the gluon two-point 
function with insertion of the CEDM operator, Fq'(p,p'), defined in Eq. (51). The third 
diagram vanishes due to the anti-symmetry of cr /liy , while the other two contribute 


a 


ryw) = V[M«“] r^(p,p') 


GG y 


2 1 e + l0g 


, m n' 1 

+ A + OlJ 


(65) 


where T^(p,p') = 4 e ,iya ^p Q p'g (see Eq. (57)). This result allows us to identify the mixing 
between the CEDM operator C and the operator O 5 = mGG. 
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Figure 4: Diagrams contributing to the gluon two-point function. The dot denotes the insertion 
of the CEDM operator. 


4.2.3 Quark-quark-gluon three-point function 

We now turn to the quark-quark-gluon three-point function with insertion of the cliromo- 


/o', - 

electric operator, T c (p,p' ,k), defined in Eq. (53). In all diagrams we chose to eliminate the 


four-momentum in favor of [k + p' — p ) M . The amputated three-point function receives 
contributions from 1PI diagrams (see Fig. [5]), non-lPI diagrams (see Fig. [2]), and quark and 
gluon wave-function renormalization. In this section we summarize our results for the 1PI 
diagrams and note that the non-lPI contributions of Fig. [2] are determined by the one-loop 


results for the quark and gluon two-point functions Y^q and Y 1 ^ presented in Eqs. (61) and 


(65), as detailed in Eq. (59). As we will discuss in Section [ 6 j we can choose a kinematic point 
and appropriate conditions so that the non-lPI diagrams are not needed to determine the 
RI-SMOM renormalization constants. 

r ^(p,p', k) can be decomposed in terms of sixteen spinor structures, and is characterized 


by sixteen scalar coefficients cp.,., 1 ^] 


~p( 3 ) _ 

1 c ~ 


Ci75 + C 2 e(e*, k,p,p ') + c 3 /*75 + c 4 ^ 75 + 75 + c 6 fl '75 


( 66 ) 


+ c 7 e M (e*, k,p)Y + c 8 e M (e*, k,p')Y + c 9 e^(e*,p,p')Y + ci 0 e M (fc,p,p , )7 / " 
+ c u a(e*, k )75 + ci 2 cr(e*,p) 7 5 + ci 3 cr(e*,p ')75 


+ c u a(k,p)'y 5 + c 15 a(k,p')^ 5 + c 16 a(p,p')j 5 


The coefficients c t are functions of the invariants p 2 , p' 2 , k 2 , q 2 , s, t, u, and e* ■ (p ±p'). The 
Ci s can be expressed in terms of triangle and bubble scalar integrals and their derivatives 
with respect to the invariants they depend on. For a generic kinematic configuration, the 
result involves logarithms and dilogarithms of ratios of invariants, and logarithms of ratios 
of the invariants to the renormalization scale fi. Working at the RI-SMOM kinematic point 
p 2 = p' 2 = k 2 = q 2 = s = u = t/ 2 = —A 2 , and in the massless limit, greatly simplifies the 
integrals, reducing them to single-scale integrals. At this point, the triangle scalar integrals 
collapse to constants, and contribute in two forms. First, triangles that are functions of three 
invariants that become equal at the renormalization point, like p 2 , k 2 and s, or p' 2 , k 2 and u, 
are proportional to the constant 

* = l (5) - H • <67) 

'Hermiticity of the operator implies constraints among the various coefficients, such as C 12 = —C 13 , which 
we have used to check our calculation. 
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The dot denotes the 



Figure 5: 1PI diagrams contributing to the quark three-point function, 
insertion of the CEDM operator. 


called Co in Ref. |59|. Here denotes the first derivative of the Digamma function. Second, 
triangles that depend on the invariants p 2 ,p /2 and t are proportional to the Catalan constant, 
which can also be expressed in terms of the first derivative of the Digamma function 

K = § (^ <I) (i) “* 2 



The only other non-rational number occurring in the result is log(2), which originates from 
the choice t = —2A 2 . 

Next, we give the UV divergent parts of the diagrams in Fig. [5j and the finite pieces of 
those Dirac structures that give non-vanishing contributions to the projections used to define 
the quark CEDM operator in the RI-SMOM scheme (see Section pi). The quark-quark-gluon 
three-point function is 


rgW,*0 


Of c 


9 <J 2a(C,fc)7 5 


Cf (£ - 2) + C A ( -j- + + log ) + h 


+ cr(e\p-p') 7s 

+ i(p + p') -£*75 


304 /I , p , _ 
^(s +l 0 S A ^] +fc2 


6C 


F - ~C A 


Clog^) + . S 


t a +... 


(69) 


where, here and later, ... denotes the contribution of the Dirac structures that are not relevant 
to defining the quark CEDM in RI-SMOM scheme. The constants k\ , /c 2 and k 3 depend on 
the definition of 75 in ci-dimension. In the ’t Hooft-Veltman scheme, they are given by 

k? v = -2 C F (2 - 0 + —+ (2 C F - C A ){ 1 - Z)K (70) 
k? V = C F {2 0 C A ~ + (2CV - C A )( 1 - 0 Q log(2) - K^j (71) 

k? v = yC F - 2C a + (-4C f + + (2 C F - CU)^y^ log(2), (72) 
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while in NDR 


kT VH = k» v -2 (c f + c a ) 

. \ i ii.' ,i~i i / 1 

rv 


NDR = k HV + - Ca 


*7 DR = - f C F + icu . 


(73) 

(74) 

(75) 


4.2.4 Quark-quark-photon three-point functions 

The quark-quark-photon three-point function with insertion of the quark CEDM gives 




with 


-C F ( \ + log ^ ^ 


+ C F cr(e*,p — p075 


+ i(p + p) -£*75 


6C f ( i + log ) + /■')" 


if )HV = C F ( -2 + It, 


kP NDR = C F ( -4 + |0 


k 


G)HV 


44 

= ( - - 40 


k 


(■ y)NDR 


= C F (8 - 4-0) . 


(76) 

(77) 

(78) 


4.3 One-loop Green’s functions with insertions of E, P, d- A and GG 


The determination of the physical block Zq of the mixing matrix in Eq. (49) requires the 


calculation of quark and/or gluon two-point functions with insertions of the operators E, 
d 2 P, (m 2 P) 12 , 3 , (md ■ A )and mGG. The renormalization of the pseudoscalar and tensor 
densities, and axial current has been studied in many papers, and the conversion between 
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The renormalization of GG 


MS-NDR and RI-SMOM to one loop was addressed in Ref. 
in MS was studied in Ref. 1601. Here we provide one-loop 1PI results for the Green’s functions 
in MS-HV and MS-NDR. 

The relevant projection of the quark-quark-photon 1PI three-point function (this is essen¬ 
tially a quark-quark function) with insertion of the quark EDM operator, evaluated at the 
symmetric point gives 


r i?’ 7) = ~ e iQA a } 


k T = C F (1 - 0 ( 2 - -0 


(1 - £)C f ( - + log ^ ) + 


+ 


(79) 


(80) 


both in HV and NDR. 
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At one loop, the 1PI quark two-point functions with insertions of the operators d 2 P and 
(t 7 i 2 P) 1 i 2 i 3 , evaluated at the symmetric point, are given by 


r<? Fi(mJp)laI = h, \ -qh\ i {VI 2 , (“} , Tr [Mh“] 1 , Tr [Mt‘\ M 


X 


47T 


kp + (3 + QCp I 7 + log Q + Cf + O 


(81) 


where kp depends on the d-dimensional definition of 75 , namely 


kp V = C f 


2(6 + 0-^^ 


„NDR _ U HV 


= kZ v -8 C f 


(82) 


The gluon and photon two-point functions with insertions of P are finite, and not needed for 
renormalization. Eq. (81) is in agreement with the result of Ref. (59], where the calculation 
was carried out in the NDR scheme. 

The 1PI quark two-point function with insertion of the operators proportional to the 
divergence of the axial current is 


r (2) _ • 
1 (md-A) 1 ,(md-A) 2 ~ 1 


Tr [Mt a ] 1 , ~ {M,t a }- -Tr [Mt a ] 1 \ -f- 


CXa 


4:/r 


k A + c F n- + io g ^ 


2 

A2 


3 

k% v = C F (£ + 4) k N A DR = C F i . (83) 

The gluon two-point function with insertion of the operator (md- A)i is finite, and we find, 


IX. = Tr M 


al ot s np 


47T 


Pa'Pg, 


(84) 


where rip = 3 is the number of flavors we are considering. The insertion of the operator 
(md ■ A ) 2 vanishes. 

Finally, the gluon and quark two-point functions with insertion of ( mGG) are given by 


p( 2 ) A 412 

(mGG) 


C A 


4Tr [Mt a ]e^p a p'^\c A 


3 + £ (1 


P 


; +l X' + fc 


a* 


6 C F ( 7 + log J 2 ) + k ° 


+ O ( M 2 ) 


fc G = -2 17-^--(3 + 0 ^ 


kc = C F (16 — 4 ijj) . 

These results determine the self-renormalization of (mGG) and its mixing with (md - A) 


(85) 

( 86 ) 

(87) 

( 88 ) 
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5 Renormalization matrix in MS scheme 


In this section we provide one-loop results for the block of the renormalization matrix given 


in Eq. (49). At various stages of the calculation we need the one-loop results for the mass, 


couplings, and held renormalization constants in general covariant gauge (recall d — 4 — 2e): 

(89) 


Z m = 1 - 3 C f 
€ 4(7T 


Z, = 1-~(C F 

€ 4(7T 

7 - 1 + 1 
7 G - — i + - — 

e 47 t 


+ Ca 


13 

~6~ 


_ 1 Oj 11 Ca — 4:T F n F 

9 e 47 r 6 


(90) 

(91) 

(92) 


where the color factors are C F , Ca, and Tp are given in Eq. (60). We will also need the 
renormalization constants for the pseudoscalar 'tp'ys'ip and tensor densities, defined by 


0{° } = Z r O r : 


Zp = ZZ. 1 = 1 + - ^ 3 C 


1 a. 


e 47 t 


Zp — 1 


1 n 
e 47 t 


(93) 

(94) 


For the mixing of dimension-5 operators, specializing Eq. (48) to the MS scheme at one loop 
one finds 


Or = (Z-‘)f Of , Z “ 5 = % - ~ % . (95) 

Note that in the above expressions a s denotes the d-dimensional renormalized coupling defined 


in Section 3.2, satisfying da s /d(\ogfi) = — 2 ea s -\-0(cc s ). So to 0{a s ) the anomalous dimension 


matrix 7 = d(log Z)/d(log/ j,) can be immediately read off Eq. (95): 7 ^ = 2 cx s /( 47 r) Zij. 
The various entries of the renormalization matrix are determined as follows: 


/n\ 

• Finiteness of the quark two-point function , gluon two-point function TjA, quark- 
quark-gluon rg } and quark-quark-photon Ep’ 7 ' ) three-point functions implies a set of 
conditions for zi n , n = 1, ...,14. Note that only the results for n = 1,..., 10 affect 
physical observables, the rest are given for completeness. 

• The operator O 2 — d 2 P renormalizes diagonally with constant Zp. 

• The quark EDM operator O^ = E renormalizes diagonally (to zeroth order in the fine 
structure constant) in the same way as the tensor quark bilinear, he., (Z ~ 1 ) 33 = Z^ 1 . 

• To zeroth order in the electromagnetic couplings, 0 4 = mFF renormalizes diagonally 
with the renormalization constant ( Z~ 1 ) 44 = Z~ x . 
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( 5 ) 

The subset of operators O 5 g 10 related to the axial anomaly renormalize, to one-loop, as 
follows 1601 (recall Z m Zp = 1): 


/ mGG \ 
(m d ■ A)i 
\ (m 2 P ) 3 ) 


MS 


/ Z^ 1 7 2 

~\%^p 0 \ 


( 

mGG 

\ 

0 

z- 1 0 



[md ■ A)i 


V 0 

0 z- 1 ) 


V 

to 

CO 

/ 


(0) 


(96) 


To explicitly check Eq. (96) at one loop use Eqs. (81), (83), (85) and (86). 

Finally, = (md ■ A )2 renormalizes as (md ■ A) 1 and = (m 2 P )\ j 2 renormalize as 
(m 2 P) 3 , thus leading to (Z-%, 88i99 = Z” 1 . 


In summary, the entries in the first row in Eq. (49) are: 

z n = 5 C F - 2 C A 
Z 12 = 0 
Z 13 = 4 Cp 
z 14 = 0 
Z15 = —2 


z 16 — Cf — -CU 

z n = 3 C F - \c A 

*18 = QCf + -C*A 

Z19 = 0 
*1,10 = 0 

£1,11 = 6CV — — Cj\ 

21,12 — —3C f + -Ca 

- ^ r 

Z 1,13 — ~^ A 

- ^ r 

2l,14 — ^C-A • 


For the remaining non-zero entries we have: 

*22 = — 3 C F 


(97a) 

(97b) 

(97c) 

(97d) 

(97e) 

(97f) 
(97g) 
(97h) 

(97i) 

(97j) 

(97k) 

(971) 

(97m) 

(97n) 

(98a) 
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^33 — Cp 

Z 44 — 3C F 

11 Ca — 4T F n F 

£55 =-^- 


+ 3 Cf 


z 56 — —QCp 

z 66 — z 77 — ^88 = z 99 = ^ 10,10 = 3 ( 7 F . 

The submatrix zu, £13 and Z 33 agrees with the original calculation of Refs. [hTI [65]. 


(98b) 

(98c) 

(98d) 

(98e) 

(98f) 


6 Definition of RI-SMOM operators and matching to MS 

A consistent phenomenological analysis of BSM-induced CP violation in hadronic systems 


requires computation of the effect of the CP-odd operators in Eq. (26) on couplings at the 
hadronic scale, such as the nucleon EDM and the T-odd nNN couplings. This is an intrin¬ 
sically non-perturbative problem. The first step in this program involves defining UV finite 
operators in a suitable renormalization scheme, whose matrix elements can be then computed 
non-perturbatively within lattice QCD. Here we will define finite operators within a class of 
regularization-independent (RI) momentum subtraction (MOM) schemes 49,59 . Next, one 


converts the matrix elements in the RI-MOM scheme to the MS scheme, commonly adopted 
to compute the Wilson coefficients and their renormalization-group evolution down to the 
hadronic scale, using continuum perturbation theory. 

In this section we address the following issues: 

1. We provide a set of regularization independent normalization conditions for the ampu¬ 
tated Green’s functions r o ( 5 ) that subtract all the UV divergences and fix the finite 


( 5 ) 

parts of the renormalization constants for the gauge-invariant CP-odd operators 0\ 10 . 

Since we will use subtraction conditions for the three-point functions at a non-symmetric 
momentum point, we call this scheme RI-SMOM, as opposed to RI-SMOM 59 . 


2. We provide the finite matching matrix that relates the RI-SMOM and MS operators to 
one-loop in perturbation theory: 


O 


RI-SMOM 


= Cij 0} 


MS 


(99) 


In practice this amounts to finding a linear combination of MS operators 0 4 MS such that 
the Green’s functions with insertions of 0^ 1-SMOM satisfy the normalization conditions 
that define the scheme (see item 1 . above). 


6.1 Defining the RI-SMOM scheme 


49,59 , with appropriate modifications related to the 


We follow the strategy outlined in Refs, 
operators we are dealing with. The content of this scheme can be summarized as follows: 
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• We require that the quark and gluon two-point functions with insertion of the quark 
CEDM operator T c\p,p') and (p, p') vanish at the symmetric kinematic point S 
defined by p 2 = p 12 = q 2 = —A 2 . 


• We require that certain projections of the three-point functions with quark CEDM in¬ 
sertion rg } and rg’ 7) take the tree-level value at the non-symmetric kinematic point S 
(involving only non-exceptional momenta) characterized by p 2 = p' 2 = k 2 = q 2 = s = 
u = t/2 = —A 2 . 


(o) 

With this choice, and by virtue of the normalization condition imposed on Y x c (p,p'), the 
non-lPI diagrams (see Eq. (59)) contributing to the three-point function with insertion 
of 0\ — C on the quark external legs vanish. In other words, the amputated Green’s 

function coincide with the 1PI Green’s functions up to a non-lPI term arising from 
operator insertion on the gluon external leg. This non-lPI term does not project on the 
spin/Lorentz structures that we use to impose the normalization conditions, so for all 
practical purposes the renormalization conditions can be imposed on the 1PI vertices. 


• We require that the gluon and quark two-point functions with insertion 0~ 6> = ( mGG) 
take their tree-level value at the symmetric point S given by p 2 = p' 2 = q 2 = —A 2 . The 
condition on the gluon two-point function involves overall factors of the quark masses. 
While one can use quark masses in any scheme, we choose to use the quark masses in the 
MS scheme. This leads to the simplest matching factors, and corresponds to imposing 
the subtraction conditions on the operator GG, ignoring the mass factors. 


( 5 ) 

• The remaining operators are related to quark bilinears: are related to the pseu- 

doscalar density, 0 \ ; is related to the tensor density, and Og_ 7 are related to the diver¬ 
gence of the axial current. We exploit this factorized structure and impose the “stan¬ 
dard” RI-SMOM conditions }59j on the quark bilinear part. The subtraction condition 

( 5 ) • • == ^ • - 

for Og_ 10 involves again overall factors of the quark masses, for which we choose the MS 

values. This is equivalent to imposing the conditions on the quark bilinears, ignoring 

the overall quark mass factors. 


Throughout, we impose the normalization conditions in the chiral limit m q —» 0. This 
is achieved as follows: (i) we expand two- and three-point Green’s functions in spin-flavor 
structures, keeping explicit powers of the quark mass, (ii) Through appropriate projections 
we then isolate the coefficients of the various spin-flavor structures, which are dehned for any 
value of the quark mass, (iii) We finally impose normalization coefficients on these coefficient 
functions in the chiral limit. This procedure dehnes a mass-independent renormalization 
scheme. 

This RI scheme, dehned in terms of gauge fixed correlation functions of quark and gluon 
states in the deep Euclidean region, serves as a useful intermediary for converting non- 
perturbative results to those required for phenomenology. In this work, we only discuss the 
matching of this RI scheme to the perturbative MS scheme in covariant gauges. To complete 
the program of connecting the MS to a lattice scheme, we also need to calculate the matching 
between lattice and this RI scheme. Among the covariant gauges, the Landau gauge is the 
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most convenient for lattice calculations. The calculation of the corresponding matrix elements 
on the lattice can be done either using lattice perturbation theory, or non-perturbatively. In 
fact, matrix elements with quark external states are used extensively nowadays for renormal¬ 
izing lattice operators |49j. However, renormalization of the CEDM operator needs extension 
of such calculations to include gluon external states. Even though gluonic correlators have 
long been studied on the lattice |66|, they are typically noisy. In addition, the matrix elements 
with two quarks and a gluon external state gives rise to “four-point” function^] and there is 
little experience with calculating these in the lattice community. 

Apart from these difficulties, however, the non-perturbative evaluation of the matrix el¬ 
ements is theoretically straightforward. The large number of off-shell operators does not 
pose a significant challenge either. In particular, since these operators explicitly involve the 
equation of motion, an n + 1-point function involving them is straightforwardly related to a 
n-point function obtained by exactly canceling an external propagator using the equation of 
motion. With such reductions, the number of correlation functions that need to be evaluated 
non-perturbatively are much fewer than the number of operators in the basis. 


6.1.1 Subtraction conditions on the quark CEDM 


We now give explicitly the fourteen conditions needed to determine Z™ SMOM 9 
the conditions on the two-point functions with external gluons and photons: 


We begin with 


e^a/JpV' 3 IT (P’P) 

e^af3P a p' 13 T c b \p,p) 


= 0 


= 0 . 


The quark-quark Green’s function has the following spin-flavor structures, 
rg } = «i 7 5 t a + a 2 cr(p,p)^5t a + «3 Mt a <J [75 + a 4 Tr [Mt a ] 

+ a 5 + a e Tr [M 2 ] 7 5 t a + a 7 Tr [Mt a ] AI 75 , 


( 100 a) 

( 100 b) 


( 101 ) 


where the cq are functions of the kinematic invariants. We impose the RI-SMOM condition 
that all the cq vanish at the symmetric kinematic point S in the chiral limit m q —)■ 0. This 
can be achieved with the following projections (traces are over color, spin, and flavor indices): 


Tr 


rS?7 5 i° 


J S 


Tr 


75 


0 

0 


ii/4- 1 


Tr 


Tg } 75^ Mt a 


Tr 


r 


Tr [Mt a 


= 0 


( 102 a) 

( 102 b) 

( 102 c) 


8 It is conventional in the lattice literature to count the point of operator insertion. _ 

9 This is in addition to the condition for eliminating possible power divergences (see Section 3.1.11. 
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/ 


M3- 1 


Tr 


rg ] 75 MH a 


Tr 

Tr 


rg } 75 t a 

rg } 75 m 


Tr [M 2 ] 
Tr [Mt a 


\ 


/ 


= 0, 


( 102 d) 


where the matrices M 2 and M 3 are given by (here Try denotes the trace over flavor indices 
only): 

Tr F [(Mt a f] (Tr F [Mt a }f 
(Tr F [Mt a }f n F {Ti F [Mt a )f 


M 2 = 


(103) 


and 


M 3 = 


/ Tr F [( M 2 t a ) 2 ] Tt f M 2 Ti> [. M 2 (t a ) 2 ] Tr F [ Mt a } Tr F [M 3 t a ] \ 

Tt f M 2 Tt f [M 2 (t a ) 2 ] (Tr F M 2 f Tr F [t a t a ] Tr F M 2 (Ti> [Mt a ]f 

\Tr F [ Mt a ] Tr F [. M 3 t a ] Tr F M 2 (Tr F [Mt a }f Tr F M 2 (Tr F [Mt a }f j 


(104) 


The above projections work for non-degenerate quark masses (m u 7 ^ rn ( i 7 ^ m s ). In the isospin 
limit rn u = rrid the matrices M 2 and M 3 become singular. In Appendix [D] we describe the 
projections needed in this case. 

To express the subtraction conditions on the quark-quark-gluon three-point function we 
restore the color and flavor indices of these objects. Recalling that T® is proportional to t a T b 
where t a is a matrix in flavor space [^] while T b is a color generator, we will use the notation 
The conditions then read (there is no summation over a and b) 


p(3) _ y p(3),ab 

c c 


1 


ie(e*,k,p,p r 

Tr 


Tr 


' { c ),ab v(p,p')t a T b 


= 2 g 


MS 


T { § ) ' ab a{kip + p')t a T 
Tr T [ c )Ab ^t a T b 


= 0 


J s 


= 0 . 


(105a) 

(105b) 

(105c) 


Note that in the first condition above, we could have used the renormalized value of the 
strong coupling constant in any renormalization scheme. The use of g MS makes the connection 
between RI-SMOM and MS schemes simpler. 

Finally, we impose the following conditions on the quark-quark-photon three-point func¬ 
tion: 


Tr 


T { c l] cr(p,p')Qt 
Tr 


= 0 


J s 


rg’ 7) 75 gt a 


= 0 . 


(106a) 

(106b) 


10 We use t° = 1/V6hx3 so that Tr F (t a t a ) = 1/2 for a = 0, 3, 8. 
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6.1.2 Subtraction conditions on the remaining operators 

We give here the subtraction conditions needed to determine the remaining entries of Zh I_SMOM . 
For the operator = ( mGG) we prescribe 


1 

6A 1 


e„v a pp a p' p vr 6 ) (p,p') 

U 5 


Tr F 


M us t a 


(107a) 


Tr 



(107b) 


(5) 

The remaining operators O 2 3 6 _ 10 are related to quark bilinears, and we wish to impose 
the “standard” 
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RI-SMOM conditions 
normalization conditions of Ref. [591 are equivalent to: 


Oo°6 _io have a simple factorized form, and the 


2 , 6—10 


A Tr 

6 q 2 


-Tr 


12g 2 


12g 2 


-Tr 


r>( 2 ) 

r o<® 75 * 




= 1 


= upTip 


J s 




o 


—Tr 
12* 

—Tr 
12 * 

—Tr 
12 * 


n (2) 

oi 5)75 


r>( 2 ) 

1 n (5)75t 


n (2) 

r '(5)75 

U 10 


= Tin 


= Tr, 


A1 ms t a 

2 \MS 


= - Ti> 
2 


= Tr, 


(7W 2 ) t a - -M MS Tr F M MS t a 
3 


(m 2 ) m V 
\m 2 )^ 

M MS t a Tr, r 


M 


(107c) 

(107d) 

(107e) 

(107f) 

(107g) 

(107h) 


The operator O 3 ; is related to the tensor density, but contains an explicit photon held 
strength. One would be tempted to impose the following condition on the quark-quark-photon 
matrix element 


12 * e(e*, k,p,p ')Tr [( Qt a ) 


Tr 


r SS °{p,p') Qt° 

U 3 


= 2e 


MS 


(108) 


J s 


which effectively fixes the projection on the structure cr(e*, k )^ 5 to its tree-level value. However, 
in terms of matrix elements of the tensor density, this prescription corresponds to 


Tr [T^{p,p')^a(p,p')\ s = 12* e^p a p' p , 
which differs from the standard one (591 

Tr [T^(p,p')v] s = 144 


(109) 


( 110 ) 
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and would lead to a finite difference in the renormalization factors [^| In our analysis we 
stick to the standard normalization condition Eq. (110). This can be obtained by imposing 


Eq. (108) while performing a finite shift Skp in the loop factor kp given in Eq. (80), namely 


5k T — Cf(1 


«I s - t* 


( 111 ) 


6.2 Matching RI-SMOM and MS operators 

We now determine the conversion matrix appearing in Eq. 

-i 


Ca = 


^■RI-SMOM\ ' 


?MS 


U 


( 112 ) 


to first order in a s . Denoting field renormalization and renormalized amputated Green func¬ 
tions of any operator O in the RI-SMOM scheme with Z q c and To, respectively, and the 
corresponding quantities in the MS scheme with Z q> c and To, the matching conditions take 
the form: 


f (2) _ 

Oi ry 




p(3) _ 

1 Oi ~ 


E c ‘i r 


( 2 ) 

Oi 


j g_ 

7j g 


E n pM 1 ' 
1 Oj 


3 

1/2 


p(3,7) _ /yy 

Oi ry 


ZqZ(j 

7 7 I /2 


E r 


(3) 

Oi 


E c ‘i r 


( 3 , 7 ) 

Oi 


(113a) 

(113b) 

(113c) 

(113d) 


When one imposes that the satisfy the RI-SMOM subtraction conditions given in subsec¬ 
tions 


6 .1.1 and 6.1.2, one obtains a system of linear equations for the CO matching factors. 


Using the explicit one-loop results of Sections |4.2| and |4.3| and the ratios of wave-function 
renormalization factors, 


Zq _ 1 . a s 

T,= ^ r,= 

Zq 47T 






(114) 


Note that in the free theory —y <j^ v ) both Eq. (1091 and Eq. ( |110[ ) hold. However, when including 


interactions a difference arises: the projection Eq. (109) selects the a component of I’^, while the projection 
Eq. (1101 picks up not only but also additional terms in E^, such as a a ^p a p'g(p^p' v — p v p'^). 
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= 1 + ? 

4n 


„ , 97 £ £ 2 \ 20 ^ 

Ca (-h —I-) — — ufTf + 

1 36 2 4 / 9 


Clf-I,- 


we solve for the C{j. 


To 0(a s ) the matching coefficients have the structure 


— r , ®s_ 

Uij — Oij -r ^ 


i I 1 

Cij + Zij log ^2 


corresponding to the RI-SMOM renormalization matrix 


7* 


RI-SMOM 


= Sii ~ An 


e + l0g X 2 ) +Q 


v 


We have given the pole terms Zij in Section [5j while the constants c 
terms of the loop factors r q c dehned above and fo,i, 2 , Au, 2 , 3 , A;) 7 3 , kc, 
in Sections |4.2| and |4.3| We fold for the first row c \ n : 


C11 = ~h - -k 2 + ^k 3 - r q 

c i 2 — 2 /o + k 2 — k 3 

r _ ill) | h~k 2 
C13 - -ffi H-2— 

C14 = 0 

C15 = -4 

_ k 2 — 2 fi 
Cu-^— 

C17 = k 2 — 2 fi 
Cis = —2/2 — k 2 — k 3 

C19 = 0 

c i,io — 0 
Ci,11 = k 2 + k 3 
k 2 + k 3 

ci,12 =-x— 


~JG 


Ci,13 — —k 2 


Cl,14 — k 


(7) 


k 3 . 


For the remaining non-zero entries of c tJ we hnd: 

C22 = ~k P - r q 


l -n F T F log^ ,(115) 

(116) 

(117) 

ij can be expressed in 
kc, and kA,p,r dehned 

(118a) 

(118b) 

(118c) 

(118d) 

(118e) 

(118f) 

(H8g) 

(118h) 

(H8i) 

(HSj) 

(118k) 

(1181) 

(118m) 

(118n) 

(118o) 
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c 3 3 = -k T - Sk T - r q (118p) 

C55 = -«g - r G (118q) 

C 56 = -kg (H8r) 

cm = C 77 = ~k A ~ r q (118s) 

Css = C 99 = Cio,io = — kp — r q . (118t) 


In Appendix [Ej we report explicit results for the matching coefficients Cy using both the 
HV and NDR prescriptions for 75 . 


7 Renormalization and the axial Ward Identities 


In the previous section we have imposed a set of subtraction conditions on CP-odd operators of 
dimension five, some of which are related to the axial current (O 6 - 7 ), the pseudoscalar density 
(^4-io)> an d GG (Of>). So far we have not discussed whether the resulting finite operators 
satisfy the non-singlet and singlet axial Ward Identities (WIs). In particular, the normalization 
conditions on the singlet A, P, and GG may be inconsistent with the singlet WIs. For the 
non-singlet case, RI-SMOM subtraction conditions have been shown to be consistent with the 


WIs [59]. 

In general one can obtain properly normalized symmetry currents through the Ward Iden¬ 
tity method, discussed in Refs. |67]68 . Moreover, in Ref. [59] the RI-SMOM conditions were 
suitably chosen so that they are consistent with the non-singlet axial WIs. Here we take a 
different point of view: we discuss how to define renormalized (singlet and non-singlet) axial 
current and pseudoscalar density operators that satisfy the axial WIs, starting from _an_ arbi¬ 
trary subtraction scheme, such as MS or the RI-SMOM scheme defined in Section 
put forward a two-step approach: 


6.1 


We 


1. Using any regulator and any subtraction scheme, define renormalized (finite) axial (A^), 
pseudoscalar (P) and GG operators. 

2. Starting from any of the above schemes, perform a finite renormalization that leads to 
operators A^, P and GG that obey properly normalized WIs. The resulting A^ is the 
“symmetry current” associated with axial transformations. We may call this new scheme 
the “WI scheme”. 


In the case of MS and the RI-SMOM scheme defined in the previous section we provide the 
explicit matching factors to the WI scheme to 0(a s ). We will also describe the procedure to 
obtain non-perturbative matching factors connecting the RI-SMOM and WI schemes. 

Our discussion is inspired by the analysis of Refs, 
theory and of Refs. 


53,69 for a dimensionally regulated 
(70,7l] for a lattice regulated theory. While we give details pertaining 
to the dimensionally regulated theory, our aim is to point out that the general features of 
the analysis are “RI” i.e., regularization independent. Therefore, we will draw parallels with 
discussions of the axial current in various lattice QCD formulations 170, 72,731 in appropriate 
places. 
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7.1 PC AC relation in terms of bare operators 

We focus on the singlet axial current for concreteness. A discussion of the non-singlet current 
in the context of dimensional regularization and minimal subtraction is presented in Ref. 1531, 
and the relevant results are a special case of the analysis presented below. In terms of suitably 
regularized operators, the PCAC relation takes the form 

d ■ A = 2(mP) + P E + X . (119) 

In dimensional regularization the bare operators take the form A M = 0(1/2) [y^, 7 5 ]0, ( mP ) = 
0 Ali 75 ' 0 , Pe = Vtsi 750 + 0 iys ipE- A" is the anomaly operator, whose tree-level insertions in 
Green’s functions vanish as one removes the regulator (d — > 4 in dimensional regularization or 
a —> 0 in the lattice theory). In the dimensionally regulated theory, with HV prescription for 
the 75 , one has 

X = ^{75, ~ , ( 120 ) 

which clearly vanishes at the classical level in d — 4 due to the anti-commutation properties 
of 75 . For d 7 ^ 4 this operator is non-vanishing and through divergent quantum corrections it 
can leave a finite remnant in Green’s functions, including anomalous terms in the axial current 
conservation equation. In NDR, X always vanishes. For this reason, NDR does not “see” the 
axial anomaly, and cannot consistently be used for the discussion of the singlet axial current. 
In the lattice theory with Wilson fermion discretization X is the variation of the Wilson term 
under axial transformation [70/721, and its properties are similar to those of X in the HV 
scheme. 

The anomalous term X can be expressed as a linear combination of other regulated opera¬ 
tors with same quantum numbers and an evanescent operator A", whose insertions in Green’s 
functions with arbitrary number of fields vanish at the quantum level as one removes the 
regulator. To perform the projection on non-evanescent operators one defines |70 

X — X + ad- A + (3 2 (mP) +7 GG , ( 121 ) 

and determines the coefficients a, (3, 7 perturbatively or non-perturbatively by requiring that 
appropriate projections of matrix elements of X in quark and gluon states (and their derivative 
with respect to the mass) vanish [^] 

0 

0 
0 

12 Insertions of X can give non-vanishing results only in Green’s functions with positive superficial degree of 
divergence. Of these, one needs to analyze only the one with two quarks and the one with two gluons, as the 
others do not provide independent information. 


(q\ x \q) 


j <75 


P-{9\x\g) 

om 


75 

(g\x\g) = 


( 122 ) 

(123) 

(124) 
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Analyzing Green’s functions of X with two quarks and with two gluons to one-loop in pertur¬ 
bation theory in the MS-HV scheme, we find 


y 2 ) 


x 


= — 4 Ci 


47T 


2075 


Ai M .75 


r 


(2 )p.v _ 


X 


= n F ^ 4 6^p a p'p , 


leading to 


a. 


ol — —4 Cp — 
47T 


an 


13 = 


7 = -n F 


(X s 

47T 


Using Eq. (121) into Eq. ( |1 19 ) one gets the final result 

(1 + a)d -A —2(1— (3)(mP) + P E - 'jGG + X , 


(125) 

(126) 

(127) 

(128) 

(129) 

(130) 


which is still expressed in terms of bare operators and couplings. The non-singlet case is now 
straightforward: one finds the same values of a and /3, and the non-singlet anomalous operator 
X a does not have a GG component. 


7.2 PCAC relation in terms of renormalized operators 

We next express the PCAC relation in terms of renormalized operators [0]j, related to bare 
operators Oj via (note that in this section we use a different notation compared to Eq. (48)) 

Oi = Zy [O], (131) 

with Zij given in an arbitrary scheme. For the operators of interest, we have the mixing 
structural 

f aa\ ( z g6 z g6m o \ ( gg \ 

(132) 


/ GG \ 


^ Z GG 

Z GG,dA 

0 ^ 


( \gg] ^ 

d ■ A 

= 

0 

Z A 

0 


[d-A] 

V ( mP ) ) 


^ 0 

0 Z m Z P y 


\ [' mP ] / 


Using Eqs. (127), (128), (129), and (130) leads to the renormalized PCAC relation 
G\{g 2 )[d ■ A] = C- 2 (g 2 ) 2[mP] + C 3 (g 2 ) [g 2 GG] + P F + X , 


with coefficients, in terms of the bare coupling g, 


Gi(g 2 ) — Za( 1 + a) + 7 Z t 


GG.dA 


(133) 


(134a) 


13 This is valid in schemes in which m Ut d,s are multiplicatively renormalized with the same constant Z m . 
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(134b) 


C 2 (g 2 ) = Z P Z m {l-P) 


Cs(g 2 ) = 


2 ^ _ 16?r2 7 7 2 7 

n F g 2 9 GG 


(134c) 


satisfying Ci i 2 , 3 ( 0 ) = 1. As a consequence of the finiteness of the EOM operator, [P E \ = P E , 
and of the independence of the operators [d ■ A], [mP] and \GG\, Ci^^{g 2 ) must be hnite. 


7.3 Finite renormalization and WIs 


Eq. (133) shows that in a given renormalization scheme (he. a choice of Za,p, Z rn , Z GG , Z gg aA 
that makes the operator insertions hnite) the renormalized quantities do not necessarily satisfy 
properly normalized (anomalous) WIs: MS-HV scheme is one example. However, given the 
scheme-dependent (g 2 ), through a hnite renormalization one can restore the WIs, as seen 


from Eq. (133). Operators in the “WI scheme” are defined by: 


Kiwi = cbo? 2 ) m 

i mP ]wi = C 2 (g 2 ) [raP} 


g 2 GG 


wi 


= Cs(9 2 ) [ 9 2 GG] 


(135) 

(136) 

(137) 


Applying the operator d/d(\ogji) to both sides of Eq. (133), using the hniteness of P E 


and the independence of the remaining operators, one obtains a set of differential equa¬ 
tions for C 1,2,3 (< 7 2 ) • The solution reveals that the coefficients C 2 , 3 (g 2 ) are such that [mP]wi 
and [g 2 GG\ wi have vanishing diagonal anomalous dimension to all orders. On the other 
hand, Ci(g 2 ) is such that [A M ] WI has an anomalous dimension starting at 0(g 4 ), related to 
the off-diagonal anomalous dimension Iqgb-a = ~ {Z~ l dZ/d{\og g)) GG d , A , namely Ja wi = 
lGGd-A ' a s/(4 7r )- The rescaled operators satisfy the properly normalized PCAC relation: 


d ■ [A] W i = 2[mP]w! + ^^]wi + + A" . 


(138) 


The coefficients needed to reach the “WI” scheme from the MS-HV scheme, to 0(g 2 ) art 14 

a 


CM. a 


C 1 = 1 - AC f -A 

47T 


C 2 = 1-8C f 

47r 


O 3 = 1 + 0(a 2 s ) . 


(139) 


In the case of the RI-SMOM scheme defined in Sec. 6.1, the perturbative values of a, f3 and 


7 are still given by Eqs. (127), (128) and (129). In HV, the conditions given in Eqs. (107c), 


(107d), and (107e), which are the equivalent to the RI-SMOM condition of Ref. 59 , give 


Of c 


Z a = 1 + -Z(AC f ) 

47r 


a. 


ZpZ m — 1 + —(8 Cf) , 

47T 


(140) 

(141) 


14 To determine C 3 we rely on the two-loop calculations of Ref. 169 . 
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where we used the value of Z m obtained in Ref. 159 

a s Cp 


— 1 


47T 




4+e-(3+e)||. 


(142) 


This leads to Ci(g 2 ) = 1 + 0(g 4 ) and C 2 (g 2 ) = 1 + 0(g 4 ), thus showing that singlet and 
non-singlet axial currents and pseudoscalar densities are already correctly normalized, up to 
corrections of 0(g 4 ). The RI-SMOM condition Eq. (107a) leads to a GG which is not correctly 
normalized. However, once a definition of Z g is given, for example by fixing the three-gluon 
or quark-gluon vertex at the symmetric point to its tree-level value |74,75] 
one to define [GGjwi- 


Eq. (134) allows 


59 


which are obtained 
1. The finite pieces in Za and ZpZ m are crucial in 


Eqs. (140) and (141) differ by a finite piece from the results in Ref. 
using NDR and found Za = ZpZ m = 
compensating the anomalous dimension of the axial current and pseudoscalar density arising 
from divergences in the MS-HV two-loop calculation, as can be explicitly verified from the 
results of Ref. 1761. For the non-singlet axial current, the cancellation is exact, and the RI- 
SMOM axial current does not have anomalous dimension at 0(a 2 ). In the singlet case, the 
0(a s ) finite piece ensures that the relation 7 a = 7 gg a-A a s/ is respected. 

While we have given explicit results in perturbation theory within the MS-HV and RI- 
SMOM scheme, the above discussion provides the steps needed to determine the coefficients 
a, f3, 7 starting from any regulator and any scheme. These, in turn, in combination with the 


renormalization factors of Eq. (132) determine the ffiiite rescaling factors 01 , 2,3 i n Eqs. (134) 
needed to obtain renormalized operators that satisfy the axial Ward identities. 


8 Relation to the AS = 1 chromomagnetic operator 


In a recent article 138 , the renormalization of the strangeness changing quark chromo-magnetic 
dipole moment (CMDM) operator has been studied. In our notation the P- and CP-even 
operator studied in |38] reads 


Ocm = g i’t AS a^ u G , 


t AS = 


/o o o\ 
0 0 1 
\0 10 / 


(143) 


Ref. 1381 studies the mixing of Ocm with lower-dimensional operators non-perturbatively on 
the lattice, and the mixing of Oqm with other dimension-5 operators in perturbation theory 
both in the lattice and in MS schemes. Clearly, a number of common issues arise in our study 
and in Ref. [38], so a closer comparison of operator basis and mixing results is desirable. 


8.1 Operator basis 

First, let us focus on the operator basis. Note that our operator basis was constructed assuming 
diagonal flavor structures, so that [Ai,t a ] = [Q,t a ] = 0 and Tr(A it a ) ^ 0. In the case of 
flavor off-diagonal generators, such as t AS , a number of new operators appears at dimension-5, 
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while all the operators involving Tr (Mt a ) vanish. In what follows we provide (i) the basis of 
dimension-5 operators mixing with the P- and CP-odd CEDM operator Oce = C defined in 
Eq. (30), with off-diagonal flavor structure t a —> t AS ; (ii) the corresponding basis for the P- 


and CP-even sector (mixing with Ocm), to be compared with Ref. |38j. 

For ease of comparison with Re 
field, i.e. 0 ( : p and ()\^ of Section 


38| we omit the operators involving the electromagnetic 
L3l With this in mind, for the AS = 1 sector we find 
ten independent CP-odd operators. In the notation of Section 3.1.3 (with t a —> t AS ), we find 


C , d 2 P, (m d ■ A) 2 , ( m 2 P)i ! 2 , Pee, d ■ A E , Ag and two new structures that vanish for diagonal 
flavor generators, namely (m 2 P) 4 = 0i7 5 [Ad, [Ad,f AS ]] 0 and (mP E ) 1 = 0 E *75 [Ad,^ 5 ] 0 — 


0 ?75 [PW, t A ‘ s '] 0 E . In order to match more closely the operator basis of Ref. 
trade the operator (md ■ A) 2 involving derivatives of the axial current in favor o 
0 E C /5 {Ad, t AS } 0 + 0 iy 5 {Ad, f A5 j 0 E , via the relation 


381 we can 

HjnP E )2 = 


2 [m d ■ A) 2 = 4 {m 2 P)i — (m 2 P) 4 + (mP E )2 ■ 


(144) 


So we end up with the dimension-5 basis in the P- and CP-odd sector reported in the left 
column of Table [5] The corresponding P- and CP-even sector operators can be obtained from 
the above ones with the substitution iy 5 —>• 1 and are given explicitly in the right column of 
Table [5] 


CP-odd 

CP-even 

Oce = ig 0(^75 G^t AS '0 

Ocm = g 0CP" G flu t AS 'ijj 

d 2 P = d 2 (0i7 5 t AS, 0) 

d 2 S = d 2 (0f AS 0) 

(m 2 P)i = \ 0?7 5 {Ad 2 , t AS } 0 

(m 2 ^)! = { 0 {Ad 2 ,t AS } 0 

(m 2 P ) 2 = Tr [Ad 2 ] 0 i 75 t A5 0 

(m 2 S )2 = Tr [Ad 2 ] 0f AS 0 

(m 2 P ) 4 = 0 ? 7 5 [Ad, [Ad,t AS ]] 0 

(m 2 ^ = 0 [Ad, [Ad,f A5 ]] 0 

Pee = #e75^ A 5 0e 

<S 0 e = 0 e^ A 5 0 e 

(mP E ) 1 = 0 i?i 75 [Ad, t AS ] 0 - 0 ? 7 5 [Ad, t AS ] 0 E 

(ms's)! = 0 E [Ad, t AS ] 0 - 0 [Ad, t AS ] 0 E 

( mP E ) 2 = 0 E i7 5 {Ad, t AS } 0 + 0i7 5 {Ad, f A5 } 0 E 

{mS E )2 = 0E {Ad, f AS ) 0 + 0 {Ad, f A5 } 0 E 

d- A e = <9 m [0£7 m 7 5 t AS 0 - 0757 / "f AS 0E] 

d -V E = i9 /i [0 E 7 M f AS '0 — 07 / d A ' s '0 E ] 

Aa = 07 5 ^t AS 0 E - 0 E $ 7 As 0 

Va = ijji(jt AS 0 E — ip E i $ t AS ijj 


Table 5: Operator basis in the CP-odd and CP-even sectors. 
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We can now compare our basis to the one in Ref. 
operators Cb . i 0 : 



which consists of ten dimension-5 


• For the gauge-invariant operators that do not vanish by the EOM, after converting 
the operators in |38] from Euclidean to Minkowski metric, we find the correspondence: 
0 1 = O cm j Oo = 2(m 2 S) 1 , 0 3 = (m 2 S ) 1 - 1/2 (m 2 S) 4 , 0 4 = -d 2 S. 

• For the operators vanishing by the EOM we find: O 5 = See, O 7 = l/2[[mSE)i — 
{itiSe) 2 ], Og = —l/2[(mSE)i + (mS#^], Og = —(d ■ Ve + Vb), O\o = Vq. 

• There is no operator in our basis corresponding to Oq in |38|. The CP-odd counterpart 
of Oq is Pee = P ) i r )b't AS ' l pEE J ripEE'ilb't AS ' l P, and it can be expressed in terms of operators 
already present in the basis, via: 

d ■ Ae = Pee + 2 Pee + (wPb )2 , (145) 

A similar linear dependence relation holds in the P- and CP-even sector. Ref. |38] finds 
at one loop that Oq is not needed to renormalize Ocm■ This is consistent with our 
finding that Oq is not linearly independent. 

• In Ref. 1381 the operator (m 2 S )2 = (m^ + m^ + mDsd is absent. This operator is allowed 
by the symmetries of the problem. In perturbation theory it can mix with Ocm starting 
at two-loop order, so its omission does not affect the results of Ref. [38]. However, the 
operator should be included in non-perturbative renormalization treatments. 

8.2 One-loop renormalization factors 

Using the CP-odd operator basis of Table [5j we have extended our analysis of the two- and 
three-point functions to include off-diagonal flavor structures and have found the mixing to 
the additional operators (in the MS scheme) 

z o CE ,(m 2 P )4 = ~ z o C e,(.™Pe)i = 74^8 ( Ca - • ( 146 ) 

Using (i) the results given in Sectio n |5| fo r the operator mixing in our original basis (extended 
to the new structures through Eq. ( ]146[ )); and (ii) the change of basis implied by Eq. ( |144 ), 
we have computed the renormalization matrix relevant to the CP-odd operators in Table [5] 

In order to compare to Ref. |38], we need the relation between the divergence structure of 
the CP-even and CP-odd sectors. At one loop we have verified that the divergences of two- 
and three-point functions with insertion of Ocm and Oce are related by a simple operation t : 

Tocm = r [r Oce\ , t : {iTs ^ 1 , t AS M n (-1 ) n t AS M n (n = 0 , 1 , 2 )} . (147) 

Similarly, the tree-level insertions of the CP-even (0 + ) and CP-odd (0_) operators appearing 
in each line of Table [ 5 ] are related by To + = t [ro_]> except for the following cases: 

1,2 n [4 1 (mP£)2,i] (148a) 
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^SeE [^p EE ] 

Fd-v E = -r [T 9 .a e ] 
r Va = -T [r^] 


(148b) 

(148c) 

(148d) 


From the renormalization matrix in the CP-odd sector and the relations (148), we have com- 

Finally, 


puted the renormalization factors in the CP-even sector, in the basis of Tabl e [5) 
converting to the basis (9^ i 0 of Ref |38] (using the relations given in Section |8T| 
our results for the renormalization coefficients to agree with Eqs. (66-75) of Ref. |3£ 


we find 


9 Conclusions 


In this work we have studied the off-shell renormalization and mixing of CP-odd dimension-5 
operators in QCD in both the MS and RI-SMOM schemes (the latter amenable to implemen¬ 
tation in lattice QCD), providing the matching matrix between operators in RI-SMOM and 
MS to 0(a s ). 

We have paid special attention to the definition of a finite quark CEDM operator in the 
RI-SMOM scheme, identifying all the needed subtractions. This is the Erst step towards 
a lattice QCD calculation of the impact of the quark CEDM on the nucleon EDM, which 
is currently afflicted by one order of magnitude uncertainty. This paper sets the stage to 
perform non-perturbative renormalization of the CEDM. The next steps in the program in¬ 
volve (i) performing exploratory computations of the needed CEDM quark and gluon Green’s 
functions on the lattice, and comparing this method to lattice perturbation theory; (ii) per¬ 
forming exploratory calculation of the CEDM insertion in the neutron state, correlated with 
the electromagnetic current or in external electric field 
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Besides inducing nucleon EDM, the quark CEDM induces T-odd P-odd pion-nucleon cou¬ 
plings that are a key input in the computation of EDMs of both light and heavy nuclei. Chiral 
symmetry implies that the T-odd pion-nucleon coupling induced by the quark CEDM can be 
extracted (up to chiral corrections) by calculating the baryon mass splittings induced by the 
quark cliromo-magnetic dipole moment (CMDM) operator 114 . In a future publication we 


will explore the non-perturbative renormalization and mixing in the flavor-diagonal CMDM 
sector and its relation to the CEDM. 

Finally, a desirable extension of this work involves studying the non-perturbative renor¬ 
malization and mixing structure of CP-odd dimension-6 operators, such as Weinberg’s opera¬ 
tor 77 and four-quark operators. 
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A CP transformation 


In this Appendix we review the definition and properties of the CP transformation. On the 
fermion fields ip, CP is defined as the linear operator CV: 

CV~ l ipCV = ip CP = *e^ 7 27oC 

= -ie^^xlY = , (149) 


where ip* = ip' T , (p is an arbitrary phase, and we are using the convention that 72 is an 
antihermitean matrix! 5 Note that the CP transformation for the t/(l)-transformed fermion 


field ip = e™ip looks like Eq. (149) with ip ip and (p —> (p + 29. 

In addition to this, the CP transformation changes all vector operators v M to v fJ in the 

metric with signature (H-), and changes every charge generator T a to {T a ) T . 

Let Ym denote a gamma structure with M Lorentz indices and 0 N denote an operator 
involving derivatives with N Lorentz indices. Then 


(xT m O n iP) cp = -{iP t O n T t m1 p X *) cp 

- OnT t m xZ {-ie-^Yox)} 

°nx 


= e iA ^ 


= e iA(p 


'ip {-^Ym^Y O n x 

= e iA ^Y c M p O NX , 


(150) 


where Y^f = (—7 2 T MX 2 Y an d we have used the hermiticity of 70, antihermiticity of 72, and 
7 q = 1. Using xj — 72* = 72; we can now write the simpler expression Y^f = —7 2 Pm 72- For 
the sixteen Clifford matrices, we then have 


1 CP = 

1 

(151a) 

7 5 CP = 

-75 

(151b) 

7? P = 

-r 

(151c) 

(7,75) CP = 

757^ = -7^75 

(151d) 

* P = 

a™ = -a*™ . 

(151e) 


For the equation of motion field ipE = ~ the transformation is 

iP% p ee (iD» CP 7m - m)ip CP 
= Xu ~ m)l2loV 

= ie‘*7 2 (-*!>* 7* ~ m)x 0 ip* 


15 Whenever we need an explicit representation for the 7 matrices, we use that one provided in Ref. 


57 
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(152) 


= - m)0* 

= ie^727o(-*-D*7 / "* - m)0* 

= ie^727o(-*-D*7 / "* - m W 
= ie l(j) 7270 [( iD^ - m)0]* 

= ie l(p ^o^* E , 

where the conjugate of D M = d fl — igA^T a is dehned as D* = d /t + igA a ^T a * to take into 
account the opposite gauge charge of the antiparticle. One way to state this result is that the 
CP phase is the same for the fields 0 and 0#, ie v = 0^,. 

Finally, note that the CP transformation on chiral helds ipL,R — (1 -F 75) / 2 "0 

CV-^lCV = ie l4, l2 {) L T 

CV-^rCV = ie^R . (153) 

B BRST symmetry and operator basis 

A given gauge invariant operator O mixes under renormalization with two classes of operators 
of same (or lower) dimension j52,;53|: (i) ghost-free gauge-invariant operators with the same 
symmetry properties of O that do not vanish by the equations of motion (EOM); (ii) “nuisance” 
operators allowed by the solution to the Ward Identities associated with the BRST symmetry. 
These include non-gauge-invariant operators. For completeness, we sketch below the procedure 
to obtain the “nuisance” operators, paraphrasing Ref. 1521. 

The gauge and fermion Lagrangian density for the SU(3)c x U(1)em group is expressed 
in terms of physical fields A“, A 7 , 0, -0, the dynamical ghosts c“, c“, c 7 , c 7 , and the non¬ 
propagating sources for BRST transformations M, M, J“, K a , J 7 , whose properties are 


summarized in Table [6j This Lagrangian is: 

C a = -07,,/-"'“' - 0 (f ■ V) 2 - ( J‘ - a„c“) !)>' ">'c h + '-,,/<•*K , ‘c i r c (154) 

- - + (d - A-<f - (7 - d,ff) 3V (155) 

+ 0 (i0 — m) 0 + M (— igc a T a — iec y ) 0 + 0 (— igc a T a — iec 7 ) M , (156) 

where 

5 *K ~ 9uA; + gf abc A\Al (157) 

F; v = d,AZ - <0A 7 (158) 

D ab c b = d,c a + gf abc Ay (159) 

= (d, - igA“T a - ieQAl) 0 . (160) 
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The action S obtained by adding to the Lagrangian density a set of infinitesimal sources <f> 
for gauge-invariant ghost-free operators O 


S — I d A xC 0 (x) + / d A x &(x)0(x) = S 0 + <f> ■ O 


(161) 


is invariant under the BRST transformations given by: 


aa; = 


dS 

aa 1 =- 6 j; sx 

(162) 

Ac a = 

8S AX 
dK a 5 

Ac 7 = 0 

(163) 

A c“ = 

^d-A a d\ 

Ac 7 = -d-A^dX 

(164) 

Axjji = 

6S AX 
dM x 5 


(165) 

A Tpi = 

8S AX 

SMi 6 1 


(166) 


with dX an anti-commuting infinitesimal parameter. This invariance leads to the Ward iden¬ 
tities for the generating functional of 1PI Green’s functions, that in particular imply the 
following identity for S = S + f d A x[ l/(2£)(<9 • A a ) 2 + l/(2£ 7 )(<9 • A 1 ) 2 ]-. 


j4 . dS dS dS dS dS dS dS dS dS dS . , 

d*x | + — —+ —+ —— I =0 . (167) 


dA^dJj dA a dJ^ dc a dK a dfcdMi dfadMi 


While S = Sq + $ ■ O satisfies the Ward identity Eq. (167), the general solution involves 
additional terms. Writing the general solution symbolically as 


S = S 0 + <S>-O + <S>-N, (168) 

and working to first order in the external sources (one operator insertion), one finds that the 
nuisance operators N must satisfy: 

W ($ • N) = 0 , (169) 


with the operator 

~ dS 0 d dS 0 d dS 0 d dS 0 d dS 0 d dS 0 d 

W = —--1---1---1---1---1-— — 

dAl dJj dJj dAl dA“ dJ“ dJ£ dA« dc a dK a dK a dc a 

dS 0 d dS 0 d dS 0 d dS 0 d 

+ WidMi + dM i lhjJ i + JfidM^ dMid^i ' 

Since WW = 0, it turns out that 

$. N = W($-F) 


(170) 


(171) 
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M 
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K 

A 

A 

4. 
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% 

Comm. 

+ 

+ 

— 

+ 

— 

— 

+ 

— 

— 

A 

Lorentz 

1 

2 

1 

2 

i 

0 

1 

2 

1 

2 

i 

0 

0 
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Color 
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3* 

8 

8 

3 

3* 

8 

8 

8 

0 

Ghost 

-1 

-1 

-1 

-2 

0 

0 

0 

1 

-1 

0 

Dim. 

5 

2 

5 

2 

3 

4 

3 

2 

3 

2 

1 

0 

2 

1 


Table 6: Properties for dynamical fields and BRST sources. The first row indicates whether 
the variable is commuting (+) or anti-commuting (—). The second and third row list the 
transformation under Lorentz and color groups. The fourth row gives the ghost number 
assignments and the fifth row lists the mass-dimension. 


where F is a set of operators with the same Lorentz property of O, same dimension, and ghost 
number —1. After acting with W one sets the sources M, M, K to zero, and to —<9 M c. 

We are now ready to classify the F operators and resulting nuisance operators N: 

• At dimension five, the only Lorentz scalars of ghost number —1 that we can write down 
are: Ax±AM , Ax±A 1 M, Mfix±Ai AA^Cxf^i Ax±$AI, M0x±4’-> McM, Mc 7 M, where 
X± = (1 i 75)/2 is a chiral projector. Acing on these structures with W produces the 
terms AeAxeA-, AeA^XeA, Ae$X±^, AAx±Ae, AA^X±Ae, A^X±Ae- 

In addition, we have the gauge-invariant ghost-free terms that are not zero by equations 
of motion in the massless limit: ^cr / “'G Ati ,x±' 0 ) A^F^XeA-i d' 2 (AXeA) , <9 /i ('0cr /J ' i/ L>*x±V ; ) ; 

d^(A^ v D v x±A)- 

• At dimension four, the only Lorentz scalars of ghost number —1 are: Mx±Vb Ax±M, 

J ail A a l0 J^A^, and Kc. The variation of these produce AeX±A, Ax±Ae, {D u G Ufl A ll + 
gAAA) - g[du,c,c]A^, (d u D u ^A'J 1 + (d^cjD^c, (D^c)c. 

The only gauge-invariant ghost-free operators not zero by equation of motion in the 
mass-less limit are: G^ U G^ V and G^G^ . 

• At dimension three and below, there are no ghost number —1 scalars, so the only oper¬ 
ators we need to consider are the gauge-invariant ghost-free operators that are not zero 
by the massless equation of motion. The only possible such terms are Ax±A- 


Selecting the T-odd and P-odd structures, including gauge-invariant ghost-free operators that 
do not vanish by the EOM, and eliminating linearly dependent operators 16 we arrive at the 


basis presented in Section 3.1 


16 We used the relation 9 Al ('0tr /il/ D v tp) = — (<9 2 + 4?n 2 ) P — d ■ Ae — 2itiPe, to eliminate one T-odd, P-odd 
structure. Moreover, there are no T-odd and P-odd operators containing the ghost fields up to dimension five. 
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C Axion Mechanism 


A very elegant way to dynamically set 9 to zero is the Peccei-Quinn (PQ) mechanism |46|, 
which predicts the existence of a new light particle, the axion 78 , 79]. We follow here the 
discussion of the PQ mechanism in the EFT framework of Ref. [801. A common feature of 
axion models is the existence of a [7pq(1) symmetry, which is spontaneously broken at high 
energy. The axion is the Goldstone boson of the symmetry, and, under [7pq(1), it changes 
by an additive constant, a —> a + c, while the SM fields are chosen to be invariant. At low 
energy, around the QCD scale, the Lagrangian includes derivative couplings of the axion to 
the quarks, which respect Upg (l). Furthermore, the symmetry is explicitly broken by the 


anomalous coupling to GG 80 , so that the quark-axion Lagrangian has the form 


e 2 a 


C = ipilpip + ]-d^ad^a + ip(C 0 + Cit 3 ) - c ail ———FF 

Ja oZTT J a 


~ [e+j)GG- M>:-„ - e-'fcWfe 


(172) 


where f a is the axion decay constant. The couplings C 0i i and c a77 are model dependent, while 
the coupling to gluons is fixed by the Ua{1 ) anomaly. 

As in Section [2J the GG term can be eliminated in favor of a complex mass term, with the 
difference that the [/^(l) rotation depends on the axion held. The rotation has the effect of 
modifying the couplings Co,i and c a77 , and, more importantly, affects the mass sector. The 
discussion of vacuum alignment of Section [2] can be immediately generalized, by replacing 0 
with 9+a/f a . In this context, vacuum alignment achieves the diagonalization of the pion-axion 
mass term. 

After imposing the vacuum alignment condition, the quark-axion Lagrangian becomes 


5 C 


~i> 


M-M - 1 


2 

mj 

2 



ip + 'ipi^s'ipm* 



(173) 


where we have kept terms quadratic in 9 + a/f a . When chiral symmetry is spontaneously 
broken, 'ip'ip acquires a non-vanishing vacuum expectation value, —(m u + m c i)('ip'ip) = 3m 2 / 2 , 
and the CP-even quark mass term in Eq. (173) generates an axion potential 

2' 


H | 0 + jJ = i(<«')Tr 


M-M "‘T 1 I » + 7 

Ja 




(■ m u + m d ) 

Vo is an even function of 6 + a/ f a , and it is minimized by 

s+V = o. 

Ja 


(m s + m d + m s - m* - [6 + — 


fa 


(174) 


(175) 
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thus canceling the CP-violating effects of the 9 term. Oscillations around the minimum de¬ 
termine the axion mass in terms of the pion mass and decay constant, and of the axion decay 
constant 

2 m\fl m u m d 

777 — --—-- 

“ P (' m u + m d y ’ 

where we neglected small corrections ~ m Uyd /rn s . 

The presence of additional, chiral symmetry breaking sources of CP violation has the effect 
of shifting the minimum of the axion potential, inducing a residual 9 term, proportional to 
the amount of CP violation. As an example, we discuss the case of CP violation from a 
quark CEDM. Performing vacuum alignment, as discussed in Section [2j induces a CP even 
axion-quark Lagrangian of the form 



5C = -|m* ) 'tpa^G^vM 1 [d C E} ^ 


+ ^ \M 1 [d C E\ — m*M 1 Tt [M 1 d C E\ } ^ + 0(9 2 ) . (177) 

When chiral symmetry is broken, the isoscalar components in Eq. (177) give a correction to 
the axion potential Eq. (174). Up to terms of [dcivi] x (9 + a/f a ) 2 that affect the value of the 
axion mass but do not change the minimum of the potential, the shifted potential reads: 

v ( s+ fyi m - (® + £) m 

- \ m * (^ + 7~) Tr [M- _1 Mce]] . (178) 

The term proportional to [dcE\, odd in 9 + a/f a , causes the potential to be minimized at a 
non-zero value of the 9 angle, 


p | <°) _ p _ r Tr r y ii _ _ {^ v gG^) 

9 + j--9 md --Tr[M [dc] j, r- ^ : 

where, by chiral symmetry, r is the same ratio defined in Section [2] 


(179) 


D Projections in the isospin limit 


We now discuss the projections needed to extract <23 ,..., 07 in Eq. (101) in the isospin limit 
rn u = rrid = m, in which the matrices M2 and M3 defined in Eq. (103) and Eq. (104) become 
singular. 

In the case a = 3, the operators O^ and 0 [g (and the structures multiplying cx 4 and a 7 ) 
vanish. To isolate 013, it is sufficient to impose 


Tr 


d 2 ) 

C 


Mt 3 


= 0 . 


J s 


(180) 
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In the isospin limit, for a = 3, a§ and are both proportional to t 3 , and cannot be disentangled 
with a flavor projection. However, the different dependence on m 2 and m 2 can be exploited, 
by imposing 


Tr 


r (2 N t 3 

\ dm 2 dm 2 ) c 5 


Tr 


d 2 

dm 2 


rg’-M 3 


0 

0 . 


(181a) 

(181b) 


The first (second) trace above isolates a 5 (c^)- 

For a = 0 and 8, M 2 is not singular even in the isospin limit, and the mixing of the CEDM 
with the divergence of the axial current is found by imposing Eq. (102c). The structures 
multiplying a 5 , a 6 and ay can be projected on the two flavor matrices t° and t 8 , so that the 
flavor projections in M3 are not independent, and the matrix is singular. Also in this case, one 
can take advantage of the different dependence on m 2 , m 2 and m s m. Defining ti = (f§ + V^o) 
and t s = (ts — ^7^0)5 «6 and 07 can be disentangled by the following projections: 


Tr 


d 2 


d 2 


dm 2 dm 2 


^ dm.dm ^ ^ 


Tr 


d 2 

dm 2 


Tr 


d 2 


ti 


dm.dm c 


= 0 


= 0 


= 0, 


where in Eq. (182a) ti (t s ) is to be used for a = 0 (a = 8). 


(182a) 

(182b) 

(182c) 


E Matching coefficients 


In this appendix we give explicit results for the matching coefficients from RI-SMOM to the 
MS-HV scheme (where i/j is defined in Eq. (67), K in Eq. (68), C A , C f and T F in Eq. (60), £ 
is the gauge parameter and n F are the number of flavors). 


C11 = 


Ca(23 + 9£) 


12 


32 C F , C A - 2 C f 

-- V> + -- 


(1 -OK 


C A ~ 2Cp 


(l + o log 2 


1 V 7> + CF (+_m + 

9 V 3 2 V 72 


-646 — 36£ + 9 £ 2 


(183a) 


C12 = (4 C F - ^(3 + 0) ^ + (C A - 2 C f )(1 — £) K + (C A - 2C F )(1 + £) log2 

+ C F ( 2 — 0 + 5 + 2£) (183b) 

ci3= (-^CV + ^(3 + 0) ip-\(C A -2C F )(l-0 K-^(C A -2C F )(l + 0 log2 
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+ c F 


(! + K> + ^ (5 - 2 « (183c) 

C 14 = 0 (183d) 

C 15 = —4 (183e) 

Ci 6 = 3c 17 (183f) 

ci 7 = ^ *l> + (C A - 2CV)(1 - 0 A' - \{C A - 2C f ){1 - 0 log 2 
o 2 

+ CV(8-0 + Cu(- 12 + 1«) (183g) 

CIS = 8CV ~^ 4(1+ ^ Ip-(C A - 2Cf)( 1 - {) A + 2 (Cg - 2 C F ) log 2 

+ C F + Pj + Cg (j - ^ (183h) 

c 19 = 0 (1831) 

c i,io — 0 (183j) 

C22 = Css = C99 = Cio,io — ~2~(3 + 0 "0 — C*f( 12 + £) (183k) 

C33 = C f (1831) 

c 55 = 2CU (l + ip + ^(-403 - 18 £ + 9 ^ 2 ) + Y n F T F (183m) 

c 56 = -4C f (4 - ip) (183n) 

c 66 = c 77 ~ ~^C F . (183o) 


The coefficients in the NDR scheme are given by (we report only the cases where c^ DR ^ 


NDR 

C 11 


NDR 

l 13 


NDR 

l 16 


NDR 

l 17 


NDR 

l 18 


r HV 

c n 


„HV 

'13 


„HV 

'16 


„HV 
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„HV 
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~3 Cf 


+ 
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Ca ~ 4 C F 

6 

C A - 4 C F 
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C 66 — C 77 


NDR 


= 0 


NDR 

c 22 
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NDR 
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= C, 


NDR 

99 


— « — "7^(3 + 0 V’ _ CV(4 + 0 
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